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FOREWORD 


This  is  an  unclassified  supplement  to  Volume  IV  of  a  series  of  reports  describing 
work  performed  by  Yale  University  under  subcontract  to  Electric  Boat  division  of 
General  Dynamics  Corporation.  Volume  IV,  Report  No.  C417-67-075,  covers  the 
period  1  July  1965  to  1  July  1966.  Electric  Boat  is  prime  contractor  of  the  SUBIC 
(Submarine  Integrated  Control)  Program  under  Office  of  Naval  Research  contract 
NOnr  2512(00).  LCDR.  E.W.  Lull,  USN,  is  Project  Officer  for  ONR;  J.W.  Herring 
is  Project  Manager  for  Electric  Boat  division  under  the  direction  of  Dr.  A.  J.  van 
Woerkom. 
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I .  Introduction 


The  standard  detector  that  one  considers  for  the  underwater  sound 
problem  approximates  a  monotonic  function  of  the  likelihood  ratio  based  on 
stationary  Gaussian  inputs  when  the  input  signal-to-noise  ratio  in  each 
channel  is  small  (locally  optimum).  Very  often  the  variances  of  the  noise 
processes  are  unknown  but  can  be  assumed  to  be  stationary  during  the 
decision  time  (quasi-stationary ).  When  this  occurs,  the  detection  threshold 
may  be  adjusted  according  to  an  estimate  of  the  noise  variance.  It  has 
been  shown  in  Report  No.  18  that  this  procedure  is  desired  for  CFAR 
detection  and  that  for  large  arrays  it  costs  very  little  in  terms  of 
detectability  in  the  presence  of  Gaussian  noise. 

It  will  now  be  shown  that  some  nonparametric  properties  are  obtained 
by  this  procedure.  That  is,  the  false -alarm  rate  can  be  fixed  for  any 
quasi-stationary  input  and  for  certain  non-Gaussian  inputs,  of  the  impulse 
variety,  the  cost  and  hence  the  miss  rate  are  reduced. 

II.  Terminology 


We  will  assume  a  threshold,  m-input,  array  detector  of  the  following 

type. 


Figure  1 


This  represents  the  locally  optimum  detector  (Report  No.  10)  for  detecting 

a  random  signal  that  is  common  to  m  channels  each  containing  stationary 

Gaussian  noise  processes  that  are  statistically  identical  but  independent. 

The  filters  are  identical  Eckart-type  filters  with  a  transfer  function 
2 

S(co)/N  (co),  where  S(<o)  is  the  signal  spectrum  and  N(w)  is  the  noise  spectrum. 
The  following  assumptions  will  be  made: 

1)  The  noise  and  signal  processes  which  are  the  components  of  the 
vector  x(t)  have  aero  mean  and  normalized  time  correlation 
functions  given  by 


P„(T) 


P,(t) 


2)  The  decision  time  T,  or  integration  time,  is  much  larger  than 
the  effective  width  of  Pn(T)»  so  that  the  test  statistic  can  be 
assumed  to  be  normally  distributed  with  negligible  error. 

When  the  variance  is  unknown,  the  detector  of  Fig.  1  is  to  be  modified 
in  the  following  way. 
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This  is  the  procedure  used  in  Report  No.  18  to  obtain  CFAR  detection  for 
Gaussian  inputs  with  unknown  variance. 

The  false-alarm  rate  (a)  of  this  detector  is  the  probability  given 
the  hypothesis  of  no  signal  present  that  S(x)  >  t(x)  .  Since  S(x)  and 
t(x)  can  be  assumed  to  be  normally  distributed  with  negligible  error,  it 
follows  that 


or,  conversely. 


r 

-h 

|s(x)-  t(x) 

jvarH|! 

3(x)-t(x)j 

4 

•  •> 
En (x +  (ji^d-a)  VarH|s(x) -t(x)j 


1/2 


(1) 


where  (jj  is  the  normalized  Gaussian  cumulative  distribution  function.  Since 

Tman2  ,  where  a R2  is  the  noise  variance,  it  follows  that 
the  multiplier  c  should  be 

[varH<js(x)-t(x)jJ  ^ 


c  ^  1  +  ^-1(l-c) 


Tma 


(2) 


As  c  ->  1  asymptotically,  it  is  sufficient  that 


c  -  1  +  $"1(1 


)  fVarH(S^  “  A  (*  ))] 
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Tmo_ 


(3) 


for  the  false-alarm  rate  to  approach  oQ  as  T  gets  large.  The  resulting 
error  in  a  is  of  the  same  order  of  magnitude  as  that  implied  in  assuming 
that  S(x)  and  A(x)  are  normally  distributed. 
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When  the  signal  is 
given  by 

2 

where  o  is  the  signal 
8 


present,  the  shift  in  the  test  statistic  S(x)  is 


variance.  The  corresponding  shift  in  t(x)  is 


Thus  for  any  m  >  2  the  adaptive  detector  is  consistent.  That  is,  the 
probability  of  detection  (p)  approaches  unity  as  T  becomes  large.  In  order 
to  calculate  the  cost  of  this  procedure  it  is  necessary  to  calculate  the 
output  signal -to-noise  ratio  of  the  detectwis,  defined  as 


and 


where 


and 


S^Rstandard 
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III.  Output  Signal -to-Noise  Ratio 


\ 


The  following  expression  can  be  evaluated  under  the  assumption  of 
stationary  noises 


VarR  S(x)  Zafaojix  [l+2dj+  amTo^h*  ,  (5) 


where 


T 

J fj  E|n12(t)  -oj*]  dt  d? 


2Tct 


(6) 


and  where  h 


•2/(1-  l)pn 


2 

(?)  d?  and  is  a  measure  of  the  "width"  of  the 


s  h' 


normalized  correlation  function  of  the  noise  processes;  and  6  •*  m  — *  — 

«  h 


2 

where  h1  is  similar  to  h  except  that  p  (?)  is  replaced  by  p_(?)  p_(?)  . 

n  ns 

The  approximation  of  Eq.  (5  )  is  a  result  of  having  dropped  a  6^  term  since 

2 

there  is  already  an  error  of  order  6  inherent  in  this  analysis,  as  a 
result  of  assumption  2 ).  Note  that  the  assumption  6  « 1  is  also 
necessary  for  the  standard  array  detector  to  approximate  the  optimum. 

This  assumption  is  valid  for  many  practical  conditions. 

When  the  signal  and  noise  processes  have  the  same  spectra,  then  b'  ■  h. 
For  Qaussian  inputs  it  is  easily  shown  that  \|r  *>  0  .  If  the  threshold  of 
the  standard  detector  were  based  on  Qaussian  inputs,  the  presence  of 
uncorrelated  noise  processes  for  which  \|r  >  0  would  increase  the  false- 
alarm  rate.  It  is  conjectured  that  the  class  of  processes  for  which  ^  >  0 
contains  those  processes  which  correspond  to  some  impulse  model.  Observe, 
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however,  that  the  dependence  of  a  on  ^  is  small  for  large  m  (m 
The  output  signal -to-noise  ratio  for  the  standard  detector  and 
is  given  by 


^standard 


»  t)  • 

for  small  6 


(7) 


For  quasi-stationary  inputs,  one  can  show  that  for  6  «  1  , 


Var^S(x)  -A(x)^  =  2m(m -lftaj4!!  1  +  2~A  6 


(8) 


It  follows  from  Eqs.  (3)  and  (8)  that  for  arbitrary  quasi-stationary  inputs. 


S^d-V 


(?) 


Therefore  the  adaptive  detector  is  asymptotically  nonparametric  in  that  a 
approaches  a  constant  (aQ)  as  T-4®  for  any  quasi-stationaiy  input.  The 
rate  of  convergence  is  limited  by  the  rate  of  convergence  of  the  test 
statistics  to  the  normal  distribution.  In  addition  one  would  expect  some 
improvement  in  detectability  for  non-Gaussian  noise  processes  for  which 

>  0  .  It  can  be  shown  that  Var^js^)- t(x)j  differs  from  VarKjs(x) -A(x)j 
only  to  order  62  provided  T  is  sufficiently  large  so  that  p  >  $0  per  cent. 
Therefore  the  output  signal-to-noise  ratio  for  the  adaptive  detector  and  for 
small  6  is 
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Vm(m-  1) 

^adaptive  ~  ,  . . . =T 

v^f 


(10) 


IV.  Cost  of  Adaptive  Procedure 

Let  us  now  equate  the  output  signal -to-noise  ratios  of  Eqs.  (7)  and 
(10)  by  considering  (m-x)  channels  for  the  standard  detector  and  m  channels 
for  the  adaptive  detector.  Thus  x  is  defined  as  the  cost  in  terms  of 
number  of  channels  for  estimating  the  noise  variance.  Solving  for  x  yields 


x  =  m  -  m-i  1 


(1-t+i)  (1-26) 


m 


m 


(n) 


For  large  m  this  can  be  approximated  by 


|  +  (t  -l)m 


2  u, 

a  h*  ^ 

SL - 1+  0 

2  .  2 
an  h 


i62 

m 


(12) 


Thus  for  veiy  small  input  signal-to-noisa  ratios,  Gaussian  inputs,  and  large 
m,  the  cost  is  approximately  ^  channel  (as  shown  in  Report  No.  18).  Let  us 
now  consider  an  example  of  a  non-Gaussian  noise  input. 

Suppose  the  input  noise  consists  of  a  process  that  has  been  obtained 
by  squaring  a  Gaussian  process  with  zero  mean  and  then  filtering  out  the 
d.c.  term.  That  is, 


I 

( 


n(t) 


Z2(t)  -CTz2 


(13) 


where  z(t)  is  Gaussian  with  zero  mean  and  variance  o„  .  It  follows  that 
n(t)  has  a  zero  mean  and  a  variance  of  2oz^. 
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From  Eq,  (6)  it  is  seen  that 
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EizU(t)aU(oi  -  l*o„W(t )z2«-)l  +  l*a  *Wz2(t )z2(C) 
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dt  dC 
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|)  PB**(e)  d£ 


-  1 


2/  (l-|)PzU(Od5 


T 

/  (1-|)PZ2(0 
•  2  +  k  ip  ■  1  ■"■  — 

I  (1-1  )Pth(Z)  d K 

o 

J.  2 

Since  pa  (£)  is  narrower  than  pa  (£)j  it  follows  that  |  >  6  .  For  example, 
if  pg(  )  is  exponential,  then  Pz^(5)  d£  ■  Pz2(25)  d£  -  ^  pa2(r])  dTl  »  and 
\jf  ■  10  ,  It  is  thus  clear  from  Eq.  (12 )  that  for  the  input  of  this  example 
x  is  negative  and  of  the  order  of  -h  channels.  That  is,  the  adaptive 
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detector  is  more  powerful  than  the  standard  detector.  This  improvement 
in  detection  is  significant  for  small  to  medium-sized  arrays.  Furthermore, 
unless  the  array  is  large  (m  »  10)  ,  there  will  be  a  significant  error 
in  the  false-alarm  rate  of  the  standard  detector. 

The  process  considered  in  this  example  has  definite  characteristics 
of  impulse  noise.  The  amplitude  density  has  a  large  peak  and  a  slowly 
(relative  to  Gaussian)  decreasing  tail.  These  characteristics  would  be 
exaggerated  if  we  raised  the  original  Gaussian  process  to  a  higher  power 
than  2  and  correspondingly  \|r  would  increase  further.  It  is  therefore 
conjectured  that  the  adapted  detector  is  more  powerful  than  the  standard 
detector  when  the  inputs  in  each  channel  are  uncorrelated  impulse  noise. 
This  effect  is  not  great,  however,  for  large  arrays  (m  »  U)  . 
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Summary 


This  report  considers  the  possibility  of  setting  a  detection 
threshold  at  each  bearing  angle  for  the  polarity  coincidence 
array.  An  adaptive  threshold  employing  the  number  of  zero 
crossings  of  the  hydrophone  inputs  is  shown  to  greatly  reduce 
the  sensitivity  of  the  false-alarm  rate  to  the  spectral 
properties  of  the  noise.  The  range  of  usefulness  of  this 
procedure  is  limited  to  the  case  when  the  self-noi3e  of  the 
hydrophones  is  such  as  to  allow  the  inputs  to  be  processed  to 
at  most  5  times  the  "nominal"  or  minimum  cut-off  frequency  of 
the  inputs.  The  analysis  assumes  that,  in  the  absence  of  a 
target,  the  hydrophone  inputs  are  all  uncorrelt  bed.  While 
this  assumption  is  unreasonable  for  submarine  arrays,  it  is 
conjectured  that  the  adaptive  procedure  can  be  modified  to 
work  in  an  isotropic  noise  field  with  unknown  space-time 
correlation. 


I .  Introduction 


This  report  deals  with  the  passive  detection  of  a  sonar  target  in  the 
presence  of  a  gaussian  noise  background  whose  spectral  properties  are 
unknown.  It  is  common  in  sonar  applications  to  display  the  output  of 
some  suitable  approximation  to  the  likelihood-ratio  detector  on  a  cathode 
ray  screen  and  have  an  observer  reach  a  decision  based  on  the  difference 
between  the  on-  and  off -target  detector  outputs.  Tiiis  procedure  is  often 
carried  out  in  spite  of  the  fact  that  the  likelihood-ratio  test  is  optimum 
only  in  the  sense  that  a  yes  or  no  decision  is  made  at  each  bearing  angle. 
There  is  of  course  no  theory  as  to  the  optimum  display,  primarily  because 
such  a  theory  would  involve  the  subjective  nature  of  the  operator. 
Psychological  studies  of  man's  decision-making  capabilities  are  just  now 
beginning.  The  decision  procedure  is,  however,  not  automated  for  a 
number  of  reasons.  The  primary  reason  is  that  the  detectors  are  based 
on  a  greatly  oversimplified  model  of  the  actual  environment  and  at  the 
present  time  do  not  have  the  capabilities  to  distinguish  between  real  and 
false  targets  that  the  operators  seem  to  have,  A  secondary  reason  is 
that  a  complete  statistical  knowledge  of  the  noise  environment  is  not 
available,  and  the  peripheral  equipment  and  computers  needed  to  measure 
_  all  the  necessary  parameters  are  too  costly  and  space-consuming. 

As  the  detectors  become  more  sophisticated,  the  desirability  of  an 
automatic  threshold  decision  made  at  each  bearing  angle  increases.  Some¬ 
times  (if  we  are  lucky)  certain  readily  measurable  properties  of  the  noise 

environment  contain  a  large  amount  of  information  in  just  the  right  form 

/ 

to  make  an  automatic  adjustment  with  minimal  peripheral  equipment.  For 
example,  given  a  number  of  simplifying  assumptions  about  the  noise 


B-l 


environment,  the  threshold  of  a  standard  array  detector  can  be  fixed  for 
non-gaussian  noise  as  well  as  for  gaussian  noise  simply  by  using  an 
estimate  of  the  noise  variance  to  set  the  threshold.  This  procedure 
assumes  that,  in  the  absence  of  a  target,  the  inputs  to  the  hydrophones 
are  all  uncorrelated  and  have  the  same  known  spectral  shape.  Of  course 
neither  of  these  assumptions  is  valid  and  a  realistic  automatic  threshold 
decision  would  have  to  take  them  into  account. 

The  main  conjecture  of  this  report  is  that  the  number  of  zero 
crossings  contains  a  great  deal  of  information  about  the  space-time 
correlation  of  the  noise  environment  and  can  sometimes  be  used  to 
automatically  adjust  the  decision  threshold.  Such  a  procedure  would 
probably  require  some  additional  processing,  such  as  pre -filtering,  bit 
hopefully  not  an  unreasonable  amount.  One  might  even  conjecture  that 
this  adaptive  procedure  could  be  nonpar ametric,  i.e,,  not  require  the 
gaussian  assumption.  However,  the  analysis  of  this  conjecture  would 
require  extensive  experimentation.  This  report  is  far  from  being  a 
complete  study  of  the  problemj  rather  it  presents  an  example  of  what 
can  be  done. 

Consider  the  use  of '  the  Polarity  Coincidence  Array  (PCA),  i.e,,  the 
Dlmus  system,  for  detecting  a  single  sonar  target  in  the  presence  of  an 
isotropic,  gaussian,  low-pass  noise  field.  Furthermore,  assume  that  in 
the  absence  of  the  target  the  hydrophone  inputs  are  all  uncar  related. 

It  will  be  shown  later  that  this  assumption,  although  reasonable  in  some 
respects,  is  for  the  purposes  of  this  stucfy  prohibitively  inaccurate  for 

# 

See  Report  No.  23  for  details. 


B-2 


\ 


feasible  submarine  arrays.  It  is  being  made,  however,  to  simplify  the 
analysis  and  to  show  a  possible  use  of  the  zero-crossing  count.  The  con¬ 
jecture  is  made  that  the  PCA  can  be  made  adaptive  relative  to  the  spatial 
correlation  by  a  similar  suitable  use  of  the  zero-crossing  count;  however, 
this  conjecture  is  not  examined.  The  method  as  well  as  the  restrictiveness 
of  such  a  procedure  are  partially  indicated  by  the  results  presented  here 
for  the  restricted  set  of  conditions. 

It  has  been  shown  in  Report  No.  22  that  for  the  conditions  just  given 
the  cost  of  clipping  is  small  when  the  inputs  are  sampled  rapidly, 
relative  to  the  cut-off  frequency  of  the  noise  spectrum.  This  cost  in 
terms  of  the  input  signal-to-noise  power  ratio  is  between  0.6  and  1.0  db, 
depending  on  the  spectral  shape.  The  cost  is  small  in  view  of  the  advan¬ 
tages  of  the  PCA,  which  are  (1)  ease  of  implementation  and  (2)  invariance 
with  respect  to  a  nonstationary  noise  power  that  varies  slowly  relative  to 
the  inverse  bandwidth  of  the  noise  spectrum.  It  was  pointed  out  that  in 
order  to  set  the  detection  threshold  it  is  necessary  to  have  accurate 
knowledge  of  the  noise  spectrum,  that  is,  of  the  cut-off  frequency  as  well 
as  the  actual  shape.  It  will  be  demonstrated  that  an  automatic  decision 
threshold  can  be  obtained  by  using  the  zero-crossing  count  in  one  or  more 
of  the  hydrophone  inputs.  This  adaptive  procedure  requires  some  elementary 
pre-filtering.  Furthermore,  this  adjustment  is  not  exact,  and  under 
certain  conditions  the  false-alarm  rate  is  appreciably  different  from  the 
designed  false-alarm  rate.  However,  in  most  cases  the  sensitivity  of  the 
false-alarm  rate  to  the  spectral  shape  is  significantly  reduced. 
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Figure  1  The  Polarity  Coincidence  Array  Detector 


II.  Terminology 

Let  us  consider  the  polarity  coincidence  array  detector  shown  in  Fig,  1, 
which  calculates  the  following  test  statistic: 

_n  r  m,  2 

3pca  '  E  I  E  iT)]  (1) 

i-1  [  3-1 

The  x^(t)  are  the  M  inputs  to  the  hydrophones,  t  is  the  sampling  interval, 

/*?  \  1  *  K  0 

sgn  (t;  =  ^  ^  <  q  ,  and  N  is  the  sample  size.  This  test  statistic  is 

compared  with  a  threshold  t,  and  the  detector  decides  that  the  signal  is 
present  if  S^ca  >  t  .  The  false-alarm  probability  a  is  the  probability 
that  Spca  >  t  given  the  hypothesis  (H)  that  the  x..(t)  are  all  independent 
gaussian  noise  processes  with  the  same  spectral  shape.  For  large  sample 
sizes  and  a  decision  time  (T)  large  compared  with  the  inverse  bandwidth 
of  the  noise  spectrum,  is  approximately  normally  distributed.  Hence 


1  -  $ 


*  - 


'VarHSp=, 


where  J  is  the  normalized  gaussian  cumulative  distribution  function,  or 

t*IH{sp4  +  3r1(1-i)VT" HW  ' 

From  Report  No.  22  jsqs,  (10 ),  (27)  and  (31  )J  it  follows  that 

eh{v}  - « !  > 
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and 


Var, 


,  4  i 

H  (  pcaj 


N  §  (H-l) 


N 

j  2 

1  +  4 

k  \  J  .  -1 

-  5f)  <;sin 

Pn(k  t) 

71 

k-l 

l 

J 

(5) 


where  Pn(T)  is  the  normalized  correlation  function  of  the  noise  inputs. 
It  is  convenient  to  define  an  equivalent  sample  size  N  as 


T  f 


N  r 

1 

Pn;icT) 

> 

n  k-l  [ 

- 

-i  > 


(6) 


where  f  is  the  sampling  rate  (f  -  -  $  N  =  T  f  ).  It  follows  from  the 
s  snrs 

previous  four  equations  that 


t 


f  T 
s 


1 


+ 


(7) 


The  equivalent  sample  size  is  shown  in  Figs.  2  and  3  as  a  function  of  the 
sampling  rate  for  a  variety  of  spectral  shapes  (See  Appendix  A).  It  is 
observed  that  N  varies  linearly  with  the  cut-off  frequency  of  the 
spectrum  (f^)  and  for  fast  sampling  varies  considerably  with  the  actual 
shape  of  the  spectrum.  Let  us  assume  that  the  second  break  in  the  spectra 
shown  in  Fig.  3  corresponds  to  the  frequency  .".ange  over  which  the  input 
data  is  to  be  processed.  It  follows  from  this  figure  that  continuous 


*The  output  signal -to-noise  ratio  of  S  a  is  proportional  to  ViT 

for  independent  samples  and  to  "\/N  1 for  dependent  samples;  hence  the 

terminology.  eq 
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operation  (f  - — >«)  is  approximately  achieved  for  sampling  rates  that  are 
s 

at  least  five  times  the  processing  range.  It  was  seen  in  Report  No.  22 
that  these  sampling  rates  are  necessary  to  achieve  the  low  cost  of  clipping. 
In  the  analysis  that  follows  it  will  be  assumed  that  the  sampling  rate  is 
sufficiently  fast  so  that  N  is  essentially  equal  to  its  asymptotic  value. 
III.  An  Adaptive  Decision  Threshold 

Let  us  now  consider  the  following  thresholds 


t'  -  f  T 


M 


1  +  l_1(l-a  ) 


K  N(0) 


(8) 


where  aQ  is  the  desired  false-alann  probability,  N(O)  is  the  total  number 
of  zero  crossings  in  a  single  channel  during  the  decision  time  (or  the 
average  over  more  than  one  channel),  and  K  is  a  suitable  constant  to  be 
determined.  The  variance  of  t'  is  quite  small  for  large  but  practical 
decision  times.  In  fact,  the  error  in  the  false-alarm  probability  caused 
by  assuming  that  t'  is  a  constant  ^N(O)  »  N(0)  j  the  average  number  of 
zero  crossingsj  vanishes  asymptotically  as  T— in  much  the  same  fashion 
as  the  error  caused  by  the  assumption  that  Spoa  is  normally  distributed. 

It  is  assumed  that  these  errors  are  small  for  practical  decision  times. 

We  will  bypass  this  detail,  however,  by  claiming  that  we  are  attempting 
to  set  the  false-alarm  probability  asymptotically  as  the  decision  time 
increases  (i.e,,  a  learning  procedure)  and  will  consider  only  the  error 
in  the  asymptotic  false-alarm  probability.  The  actual  asymptotic  false- 
alarm  probability  (o)  is  obtained  by  equating  t  of  Eq.  (7  )  to  t 1  of  Eq.  (8) 
and  replacing  N(0)  by  N(o),  resulting  in 
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(9) 


a  =»  1  -  Q  Jf-1(l-a  )  n/  ■  "e3" 

0  VKN(O) 

We  will  now  discuss  the  reasoning  behind  using  t'  as  a  detection  threshold. 
The  well-known  expression  for  the  average  number  of  zero  crossings  is  [l] 

•  n1/2 


N(0)  =  2T 


/ 


r  s(f )  df 


I  S(f)  df 
0 


(10) 


where  S(f )  is  the  spectral  shape  of  the  noise  inputs.  Let  us  assume  that 
S(f )  can  be  written  in  polynomial  form  and  that  there  exists  a  first  break 

f 

point  in  the  spectrum  denoted  by  f. .  If  we  set  y-  -  x  ,  we  can  write 

1  h 


1/2 


N(0)  »  2Tf^ 


(11) 


This  expression  can  be  readily  evaluated  for  a  variety  of  low-pass  spectral 

shapes  (See  Appendix  B).  Recall  [l]  that  for  a  "single-pole"  spectrum, 

N(0)  ■  «  ,  and  our  system  obviously  fails  (a-+$0  per  cent).  We  will 

assume  that  there  is  a  processing  filter  somewhere  (maybe  the  hydrophones 

themselves  )  that  prevents  this  possibility. 

Figure  U  is  a  plot  of  N  (0 )  versus  N  (large  sampling  rates  )  for  a 

eq 

variety  of  spectral  shapes  that  can  be  obtained  by  passing  white  noise 
through  a  low-pass  filter  having  no  zeros  and  only  real  poles.  For  these 
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spectra  the  following  conclusions  can  be  reached.  For  spectra  that  fall 


off  sharply 


r 

- 

S(f)  - 

1 

i  n>3 

1  + 

f 

,  there  is  very  nearly  a  linear 


relationship  between  N  and  N  (6).  For  spectra  that  fall  off  slowly  but 

eq 

N 

-22-  curve 


are  then  cut  off  sharply  at  some  frequency  f ^  >  f^  ,  the 


nTo) 


has  the  same  linear  relationship  for  fg  <  2f^  and  falls  off  somewhat 
for  larger  fg. 

It  follows  that  we  should  pre-filter  the  noise  inputs  with  a  filter 
having  the  transfer  function 


H(» 


H 


1  + 


“2 


3 


(12) 


where  is  equal  to  the  processing  frequency  in  radians.  The  constant  K 
needed  for  the  decision  threshold  [Eq.  (8)J  is  seen  from  Fig.  U  to  be 
2.7UO  .  If  the  noise  spectrum  begins  to  fall  off  slowly  at  some  frequency 
f^  considerably  less  than  fg»  then  there  will  be  an  error  in  the  asymptotic 
false-alarm  probability  given  by  Eq.  (9 ).  This  error  will  be  examined  in 
more  detail  in  the  next  section.  In  addition  to  the  spectra  indicated  in 
Fig.  U,  the  presence  of  zeros  in  the  spectrum  (filter  transfer  function) 
will  also  be  analyzed  (spectra  restricted  to  be  monotonically  decreasing). 

It  is  conjectured  that  the  suboptimum  PCA  discussed  in  Report  No.  22 
will  behave  similarly  with  respect  to  the  adaptive  threshold  and  in  fact 
be  somewhat  better  for  large  arrays.  This  conjecture  is  based  on  the 


fact  that  N  of  the  suooptimum  device  has  a  larger  relative  increase  thaa 
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the  optimum  PCA  as  the  noise  inputs,  prior  to  processing,  vary  from  white 
to  one  having  a  single-tuned  spectrum.  It  follows  that  the  contours  of 
Fig.  U  will  probably  be  even  closer  to  a  straight  line  for  the  suboptimum 
PCA. 

IV.  Error  in  Asymptotic  False-Alarm  Probability 

The  actual  asymptotic  false-alarm  probability  can  be  determined  from 
the  designed  false-alarm  probability  aQ  and  Eq,  (9)  using  Fig.  U  (or  the 
data  in  the  appendices ).  It  will  be  assumed  that  the  input  data  are 
processed  according  to  the  previous  section.  The  false-alarm  rate  is 
shown  in  Fig.  5  for  a  variety  of  spectral  shapes  and  a  designed  false-alarm 
rate  of  1  per  cent.  The  abscissa  is  f^/fg  *  where  f^  is  the  cut-off 
frequency  of  the  noise  spectrum  and  fg  is  the  processing  frequency.  As 
indicated  previously,  the  only  difficulty  arises  when  the  noise  spectrum 
falls  off  slowly  and  is  significantly  reduced  by  the  time  the  processing 
frequency  is  reached.  The  least  favorable  case  ("single -pole"  spectrum) 
is  shown  again  in  Fig.  6  for  different  designed  false-alarm  rates. 

Thus  we  see  that  some  knowledge  of  the  spectral  shape  is  necessary 
for  an  accurate  setting  of  the  false-alarm  probability.  Suppose  for 
instance  that  we  know  that  the  spectral  shape  is  always  flat 

below  some  frequency  f^.  The  maximum  error  in  a 

caused  by  processing  the  inputs  out  to  fg  =  9f^  can  be  obtained  from 
Fig.  6  for  the  value  f^/fg  =  0. 14  .  This  is  quite  a  reasonable  maximum 
error.  Sometimes  one  has  some  other  specific  knowledge  of  the  noise 
spectrum  that  enables  one  to  make  an  alternate  setting  of  the  threshold. 
Consider  the  following  example. 


[s(o)-s(f)] 

S(0) 


<  0.2 
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Let  us  suppose  that  the  spectrum  commonly  appears  to  be  "single-tuned" 
with  the  cut-off  frequency  in  the  vicinity  of  1250  cps.  Furthermore,  let 
us  assume  that  the  self-noise  of  the  hydrophones  is  such  that  processing 
out  to  5  kc  seems  reasonable.  Finally  assume  that  the  noise  spectrum 
never  falls  off  faster  than  12  db/octave.  We  wish  to  set  the  false-alarm 
probability  automatically  in  such  a  way  that  gross  errors  do  not  occur 
for  noise  spectra  other  than  the  nominal  one.  We  therefore  process  out 
to  5  kc  and  set  the  threshold  to  give  no  error  for  the  nominal  spectrum 
(K  -  2.1U& )  .  The  resulting  actual  false-alarm  probability  is  shown  in 
Fig.  7  for  a  designed  false-alarm  rate  of  1  per  cent.  Included  in  this 
figure  is  the  actual  false-alarm  rate  for  the  non-adaptive  detector  whose 
threshold  is  based  on  the  nominal  spectrum.  We  can  reach  the  following 
conclusions . 

For  the  non-adaptive  detector,  the  uncertainty  in  the  spectral  shape 
causes  a  possible  variation  in  a  of  5  db.  A  similar  variation  is  caused 
by  a  20  per  cent  uncertainty  in  the  cut-off  frequency  (1  kc  <  f ^  <  1.5  kc)  . 
For  these  uncertainties  in  the  spectral  shape,  the  threshold  of  the  non- 
adaptive  detector  can  be  set  to  guarantee  a  maximum  error  in  a  that  is 
less  than  +  5  db.  However,  the  threshold  of  the  adaptive  detector  can 
be  set  to  reduce  the  maximum  error  in  c  to  i  0,7  db  for  the  same 
uncertainty  in  the  spectrum.  Figure  8  is  a  similar  plot  with  the  only 
difference  being  that  the  inputs  are  processed  only  out  to  the  nominal 
cut-off  frequency.  With  this  type  of  processing  the  false-alarm  rate 
becomes  relatively  insensitive  to  the  spectral  shape.  As  a  result  the 
non-adaptive  detector  threshold  can  be  adjusted  to  guarantee  a  maximum 
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error  in  a  of  only  -  1.5  db  for  a  20  per  cent  uncertainty  in  the  cut-off 
frequency.  On  the  other  hand,  the  adaptive  detector  introduces  virtually 
no  error  for  this  type  of  processing. 

V.  The  Effect  of  White  Self-Noise 

Finally  we  will  consider,  in  a  limited  way,  the  effect  of  allowing 
a  zero  in  the  spectrum.  More  specifically  we  will  consider  a  spectrum 
of  the  following  type:  r 


S(f)  - 


(13) 


where  f^/f^  <  1  and  f^/^  <  1  .  Hence  f^  is  still  the  first  break 
point  in  the  spectrum.  This  spectrum  can  alternatively  be  written  as 


S(f) 


-51 

p +  q 

- 

j 

. ,!.• 

1  + 

f 

1  + 

f  h 

\T1 

0  J 

klj 

(1U) 


where  q  is  a  constant  equal  to 


1  - 


3'  J 


.  Thus,  if  the  second  term 


is  regarded  as  the  processing  filter,  the  zero  of  the  spectrum  corresponds 
to  a  minimum  level  in  the  noise  power  spectrum.  Small  values  of  q  would 
always  be  present  as  a  result  of  the  self-no.i3e  generated  in  the  hydrophones. 
In  fact,  in  the  absence  of  specific  knowledge  of  the  signal  spectrum,  the 
inputs  are  normally  processed  out  to  some  frequency  determined  by  this 
self -noise  (f g  <  f^)  in  the  hope  that  the  signal  spectrum  is  similar  V 
the  incoming  noise  spectrum. 
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Let  us  now  assume  that  the  adaptive  threshold  is  adjusted  under  the 
assumption  that  the  noise  spectrum  is  white  (prior  to  the  processing 
filter)  or  K  -  2.3U55  .  The  actual  false-alarm  rate  (oq  ■  1  per  cent) 
is  shown  in  Fig.  9  for  a  number  of  locations  of  the  pole  where  the 
abscissa  represents  the  location  of  the  zero.  For  small  values  of  self¬ 
noise  (f^/f^  small)  ,  the  spectrum  falls  off  at  a  rate  slightly  less  than 
6  db/octave,  and  hence  the  error  in  a  increases.  Since  the  threshold  was 
designed  to  give  no  error  for  white  noise ,  the  error  starts  to  decrease 
when  the  self-noise  becomes  dominant  in  the  vicinity  of  the  processing 
frequency.  Both  effects  will  be  far  less  noticeable  for  speotra  that  fall 
off  faster  than  6  db/octave.  No  analysis  has  been  made  for  higher-order 
zeros  or  for  a  simple  zero  where  f^/f^  >  1  (Spectrum  increasing  in 
some  frequency  range). 

Let  us  re-examine  our  last  example  (nominal  cut-off  frequency  of 
1200  cps ).  Let  us  now  assume  a  self-noise  -20  db  that  of  the  low-frequency 
noise.  Since  the  spectrum  can  fall  off  as  fast  as  12  db/octave,  it  seems 
reasonable  to  process  out  to  5  kc.  Figure  9  tells  us  (f^/fg  -  0.25  > 
f^/f^  ■  ,1 )  that  for  the  assumed  spectral  shape  the  false-alarm  rate  is 
increased  by  less  than  1  db.  For  sharper  spectra  this  increase  is  smaller. 
It  follows  that  the  previous  decision  to  set  the  threshold  for  the  nominal 
spectrum  (no  white  noise)  is  about  right.  This  will  insure  a  maximum  error 
in  the  false-alarm  rate  of  about  +  1  db  within  the  tolerance  limits. 

It  appears  that  whether  we  are  setting  the  threshold  on  the  basis  of 
some  lowest  cut-off  frequency  or  some  nominal  cut-off  frequency,  the 
adaptive  procedure  works  well  for  processing  frequencies  as  large  as  U 
times  the  nominal  (lowest)  cut-off  frequency.  For  reasonable  amounts  of 
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self-noise,  there  does  not  appear  to  be  much  advantage  anyway  In  processing 
to  seme  higher  frequency.  In  fact  it  is  the  presence  of  self-noise  that 
makes  the  adaptive  system  work.  For  with  extremely  small  amounts  of  self¬ 
noise,  one  should  process  the  inputs  to  some  relatively  high  frequency. 

If  one  processes  to  some  frequency  fg  >  10  f^  ,  then  the  adaptive 
detector  is  no  better  than  the  non-adaptive  detector. 

VI.  On  the  Assumption  of  Uncorrelated  Hydrophone  Inputs 

Let  us  assume  a  linear  array  of  hydrophones.  The  mean  of  the  PCA 
detector  output  is  given  by 


(15) 


where  p  (d)  is  the  normalized  spatial  correlation  between  two  hydrophones 
spaced  d  feet  apart.  It  follows  that  the  threshold  [Eq.  (7 )]  should  be 
modified  to  be 


t  -  f  T 
8 


M 


<[i+r]  +  f^Cl-tt) 


(16) 


M 

where  R  -  ^  (1  -  j-j)  sin"1  j^pg(kd)J  ,  and  is  somewhat  smaller  than 

k-l 

N  as  a  result  of  the  spatial  correlation.  We  wish  to  indicate  that  the 
eq 

assumption  of  statistical  independence  (R  -  0)  is  prohibitively  inaccurate 
for  feasible  submarine  arrays,  at  least  as  far  as  setting  a  detection 
threshold  is  concerned.  To  do  this  we  assume  that  can  be  accurately 
determined  (say,  by  the  zero-crossing  count ).  It  can  easily  be  seen  that 
the  error  in  the  false-alarm  rate,  due  to  R  /  0  ,  can  be  calculated  from 
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I  in  1 

r1(l-a)-r1(l-aft)-1/ - r  (17) 

0  V2(l-§) 

Let  us  now  place  the  hydrophones  sufficiently  far  apart  so  that  the 

error  in  a  caused  by  the  spatial  correlation  is  "only"  +3  db  for  a  designed 

false-alarm  rate  (aQ)  of  1  per  cent,  jin  reality  this  error  is  very 

insensitive  to  a  and  varies  from  2  db  to  i*  db  as  o  varies  from  10  per  cent 
o  o 

to  .01  per  centfj  We  will  call  this  spacing  the  minimum  spacing  for 
"uncorrelated"  inputs.  This  spacing  is  achieved  when 


<  0.25 


(18) 


It  has  been  shown  in  Report  No.  1  that  for  a  three-dimensional  isotropic 
noise  field  and  a  single  pole -spectrum, 


Ps(kd) 


< 


c 

2nfxkd 


where  c  is  the  speed  of  sound  in  ft/sec.  This  upper  bound  is  very  nearly 

achieved  for  small  values  of  p  (kd)  and  for  arrays  that  are  not  steered  in 

s 

the  end-fire  condition.  For  these  conditions 


M 


*  2  .cir 
7  * i  *  ^ 


(i  _  *)  i 


k“l 


_1_  c  i  \k 
2n2  *1  d  12 


+ 


*  M  iil 
2?  d 


(19) 
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where  H  is  very  nearly  equal  to  M  for  small  values  of  M  (M  <  6)  and  is  less 
than  M  for  large  arrays.  It  follows  that  the  "minimum"  spacing  is  given  by 


The  equivalent  number  of  uncorrelated  samples  is  normally  quite  large, 
resulting  in  a  large  d.  For  lightly  correlated  inputs  with  single-tuned 
spectra,  N'  is  of  the  order  of  10 Tf.  (see  Fig,  2  or  3).  If  we  now 
assume  that  f^  <  5  kc  ,  a  lower  bound  for  the  minimum  spacing  becomes 

d  >  15  — ...J - YT1  .  (21) 

p 

The  function  — .......  varies  frcm  2  to  10  as  M  varies  from  2  to  50. 

V2(1-a> 

The  minimum  corresponds  to  M  -  2  and  hence  d  >  90  Vt"1  ft.  Thus  for 
wide  spectra  (f^  ■  5  kc )  and  short  decision  times  (T  1  sec ),  an  array 
with  2  hydrophones  spaced  90  ft  apart  has  crudely  independent  inputs.  On 
the  other  hand,  for  a  six-element  array,  d  >  200 ft  »  and  hence  one 
needs  a  1000-ft  array  even  for  short  deoision  times.  We  conclude  that  the 
spatial  correlation  of  the  noise  field  must  be  considered  for  practical 
submarine  arrays  in  order  to  set  a  decision  threshold  independently  at 
each  bearing  angle. 

It  is  speculated,  however,  that  a  similar  procedure  can  now  be 
carried  out  to  adapt  for  the  spatial  correlation  as  well  as  the  time 
correlation.  For  instance,  it  would  seem  reasonable  to  try  an  adaptive 
threshold  of  the  following  type: 
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*"  -  fgT^ 


(22) 


•s 


Tf(M) 

1  +  - ± - +  1(1-  a) 


f,(M) 


K,  N(0) 


■»/: 


where  f^(  )  and  fg(  )  are  suitable  functions  and  and  Kg  are  suitable 
constants.  It  should  be  pointed  out  that  the  variance  of  t"  is  no  longer 
negligible  and  will  have  to  be  considered.  However,  for  large  arrays,  if 
N(0)  represents  the  average  number  of  zero  crossings  over  all  the  inputs, 
the  cost  introduced  by  the  variance  of  t"  should  be  snail.  A  check  of 
the  validity  of  this  modified  adaptive  threshold  would  be  quite  cumbersome 
and  has  not  been  carried  out  to  date. 
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Appendix  A  Calculation  of  N 


£3 


was  defined  according  to  Eq.  (6)  as 


T  f 


eq 


N 

a  - 1)  ■ 

sln^UkT)! 

n  k-l 

L  JJ 

where  Pn(T)  is  the  normalized  correlation  function  of  the  ncise  inputs. 

Given  p  (T),  N  can  therefore  be  calculated  directly  on  a  computer.  The 
correlation  function  is  obtained  via  the  Wiener-Khintchin  theorem  from 
the  power  spectrum: 

p(T)  =  j S(f)  cos  2nf  df  ,  (A-l) 

0 


where  p  (t)  = 

n  P(0) 


Assuming  a  first  break  point  in  the  spectrum 
a  f 


occurring  at  f.  and  setting  x  -  -jr-  ,  one  obtains 

*1 


p(kr)  »  S(x)  cos  xy  dx  , 


(A-2) 


a  2nkfl 

where  y  -  — — —  .  We  will  consider  three  classes  of  power  spectra. 
fs 

Case  1.  Consider  spectra  of  the  form 


S(x) 


1  +  x 


3H 


(A-3) 


where  n  is  any  positive  integer.  Equation  (A-3)  can  be  evaluated  directly 
[2],  resulting  in 
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p(kT) 


y  cos 


lA-U) 


where 


'('i-U  Z>ITrJ’1 


The  limiting  case  (n— »<*>)  is  of  course 


p(kT)_+fili2J[  . 


(A-5) 


The  calculation  of  Eq.  (6)  is  then  straightforward  and  is  plotted  in  Fig.  2 

for  n  ■  1,2,3.  and  <*■  versus  the  normalized  sampling  rate  f  /f,  .  Observe 

that  for  f ^f^  =  100  ,  N  has  essentially  reached  its  asymptotic  value. 

We  will  therefore  use  this  calculation  of  N  in  order  to  set  the  threshold. 

eq 

Below  is  a  table  of  these  results. 


N^T^  11.530  U.691  3.883  3.163 


Case  2,  Consider  spectra  of  the  form 


/,  2  w.  2n  2n 

(1  +  x  )(l  +  p  x  ) 


(A-6) 


where  p  ■  f^/fg  and  n  >  1  .  This  can  be  factored  into  the  form 


B-27 


’i-<-iry"1  [1+ 


,  2n(  2n-2  2n-U  2n-6  ,  -  \n+l| 

1  p  lx  -x  x  -  (-1)  I 


,  ,  2n  2n 
1+  p  x 


(A-7 ) 


This  can  now  be  transformed  term  by  term  with  the  correlation  function 

corresponding  to  S'  (x)  = - 5n”Sn  ^or  ^  <  m  <  2n-l  b©ing  0  for  m  even 

1  +  p  x  n 

and 


p'(kT) 


n  y  .  f  (2k-i )  1 

tfiy’e'p  1!"  "jslnjaL^„+zCos%^„ 

1  2n  m  2n  P  2n 


(A-8) 


for  m  odd.  Note  that  for  n  odd  and  greater  than  one,  this  procedure  leads 
to  an  indeterminant  form  for  the  value  p  -  1  .  For  other  values  of  p 
there  is  no  difficulty. 

For  example,  let  us  assume  that  n  =  2  .  Then  it  is  readily 
determined  from  Eqs .  (A-7  )  and  (A-8 )  that 


p<kT>  ■  fi  5  rrHr  e'y*  °  pV?  p3  sl" ! + 


and  hence 


%(kT)  '  3  1 

1*  — (p3-p)J 


-y  .  '  „Vj 


e  +  e 


*e  P  2  p3elnF wt'P0MFw?  ‘ 


(A-9) 


Thus  for  a  given  p  Eq.  (6)  can  be  calculated  in  a  straightforward  manner. 

As  before,  the  value  for  N  at  f  /f.  »  100  will  be  assumed  to  be  the 

eq  s  1 

asymptotic  value  for  N  .  We  will  now  tabulate  the  results  (N^Tf^ )  for 
a  variety  of  values  of  n  and  p. 
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\p 
n  \ 

.1 

.25 

.5 

— 

1 

5 

10 

1 

9.133 

7.21*65 

5.630 

b.103 

— 

— 

2 

9.358 

7.198 

5.230 

3.1*06 

U.U8 

ii.62 

3 

9.ii01 

7.210 

5.132 

— 

— 

— 

Case  3.  Consider  spectra  of  the  form 

S(x)  - - -  (A-10) 

(1+  x  )(1+  p  x  ) 

This  can  alternatively  be  written  (provided  q’  <  1  ,  which  guarantees  that 
f^  is  still  the  first  break  point  in  the  spectrum)  as 


S(x)  c< 


(1  +x2)(l+  p2nx2n)  (l  +  p2nx2n) 


(A -11 ) 


It  is  clear  that  the  correlation  function  can  now  be  obtained  by  a  linear 
superposition  of  case  1  and  case  2.  Hence  for  n  =  2  , 


p(kT)  -  f  j  -±-r 
1+P 


. 

— 

\ 

T  *  -Z-l  ■ 

"y  pVi1 

eft  4>  p  “  T 

p^ +  q(p^+ 1)  3in(n  y  \  D  co3 

S  D  “  C 

p  l1  V?p  1 

1*  v?p| 

/ 

and  hence. 


PnOtO 


“y .  pV¥ 


1+^p3-p+aiE^] 


<e  +  e 


1 3in 


l* 


-  P  cos  hj 


?  +  -Z- 


V^P  I 


_Z_ 

VFPi 

(A-12) 


We  will  now  tabulate  the  results  (N  /Tf-  )  for  the  case  n  =  2  . 

eq  1 
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\q 

p  \ 

0 

.001 

.003 

.01 

.05 

,10 

.5 

1,778 

GO 

n 

9.358 

9.1*85 

9.750 

10.705 

15.512 

20.1*36 

36.21*8 

— 

1*6.910 

.25 

7.198 

7. 2U0 

— 

7.350 

8.565 

9.679 

114.177 

— 

18.761* 

.5 

5.230 

— 

— 

5.308 

— 

5,888 

7.279 

8.1*51 

9.382 

Appendix  B  Calculation  of  N(0) 

The  average  number  of  zero  crossings  can  be  calculated  directly  from 
the  spectral  shape  from  Eq.  (11). 

-1/2 

r  cd  i  * 


N(0)  =■ 


(B-l) 


This  expression  can  be  evaluated  readily  from  the  results  of  the  previous 
appendix  simply  by  recognizing  that 


S(x)  dx  -  p(0) 


and 


(B— 2 ) 


S(x)  dx 


p»(0) 


> 


0 


where  p*(0)  is  given  by  Eq.  (A-8),  where  m  is  suitably  defined.  Let  us 
explicitly  evaluate  N(0)  for  the  same  classes  of  power  spectra  considered 
before. 

Case  1.  For  power  spectra  of  the  form 

S(x) - —  ,  (B— 3 ) 

l  +  x2n 

wo  obtain  from  Eqs.  (A— U )  and  (A-8)  that 

k^l 
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and 


‘"<°>-fisZ  sl”Iiri23" 

k«l 


where  the  latter  result  is  valid  only  for  n  >  2  ^Recall  N(0)  -  » 
"single-pole"  spectrum.!  .  Using  the  identity 


n 

£  sin(2k  -1)  0 
k-l 


sin2  n9 
sin  0 


) 


(B-U) 


for  a 


one  readily  evaluates  Eq,  (B-l)  as 


N(O) 


2T  f. 


Sin£ 

■*¥ 


1/2 


» 


for  n  >  2  .  The  results  are  tabulated  below. 


n 

2 

3 

12 

flO 

nToJ/t^ 

2 

l.iau 

1,216 

1.155 

(B-S^ 


Case  2.  For  power  spectra  of  the  form 


S(x) 


_ 1 _ 

d+x^d+p^x211) 


(B-6) 


p(O)  is  directly  determined  from  Appendix  A;  however,  a  different  factoring 
procedure  is  needed  to  evaluate  p'(0). 


x2S(x) 


i*(-i  rVn 


i  ,  i  *  A2*-2  -  * . . . .  (.1  r  w 

.2 


1  +  x 


-  2n  2n 
1+  p  x 
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One  can  now  use  Eq.  (A-8)  to  determine  p '(°)-  For  example,  let  us  assume 


that  n  -  2  .  Then 


and  from  Eq.  (A-9), 


p»(0)  -  f  |-1  +  p^|-  +  (  > 


p(o)  *»  |  ^Tiji  +  TT  "  p  “T"  r  * 


(B-8) 


(B-9) 


Hence  from  Eq.  (B-l), 


N(0)  «  2TJT.  -v  ItJB.-Z.l&'p. 

VpU-p2+V7p 


(B-10) 


These  results 


n  and  p. 


NCoj/T^ 


will  now  be  tabulated  for  a  variety  of  values  of 
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(B-12  ) 


p( 0)  -  f  5— L-di  +  ^ 
i  +  p 


p3-p+j|  (p^+D 


Fr cm  Eqs.  (B— 8 )  and  (A-8)  we  readily  obtain 

pf(o)  -  5 


n  1  J  ,  .  V? 


W 


-i  +^r  (p+p}r +  fi?^ 


2.V? 


-I)  + 


,3+p+a  (p^* + 1 ) 


Hence  from  (B-l ), 


N(0) 


/  P3+P+f  (p^+l)]  -  P2 
-  2  T  f ,  - i-i - 2 - i - 


\F  2 

T  p 


j3-p+3  (p^  +  l)  +  p2 


We  will  now  tabulate  the  results 


NfO/T^ 


for  the  case  n  »  2  , 
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Summary 


The  ambient  noise  model  proposed  by  Talham  (JASA  36,  151*1, 
196U)  is  used  to  estimate  the  vertical  directionality  of  ambient 
noise  for  a  variety  of  velocity  profiles  and  hydrophones  located 
100-500  feet  below  the  ocean  surface.  Except  when  the 
hydrophone  is  located  in  a  region  of  negative  velocity  gradient, 
the  model  indicates  that  much  of  the  ambient  noise  arrives  from 
a  near-horizontal  direction.  The  calculated  difference  between 
horizontal  and  vertical  intensities  exceeds  20  db  in  a  number 
of  cases.  When  the  hydrophone  is  located  in  a  region  of 
negative  velocity  gradient,  there  is  a  null  in  the  noise 
directivity  pattern  in  the  immediate  neighborhood  of  the 
horizontal.  Aside  from  this  obvious  effect  of  ray  geometry, 
the  pattern  is  not  markedly  different  from  that  of  a  hydrophone 
located  in  a  region  of  positive  velocity  gradient. 
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I.  Introduction 


In  underwater  sound  detection  systems,  ambient  noise  plays  a  major 

role  in  limiting  target  detectability.  It  has  been  suggested  by  several 

researchers  that  a  large  contribution  to  the  total  ambient  noise  is 

generated  at  the  surface,  both  by  wave  splashing  and  by  wind  blowing 

across  wave  tops.  In  particular,  R.  J.  Talham^  has  used  a  model  in 

which  the  ocean  surface  is  considered  planar  and  covered  by  a  uniform 

distribution  of  independent  noise  sources.  The  purpose  of  that  paper  was 

to  determine  the  vertical  distribution  of  ambient  noise  as  received  by  a 

hydrophone  located  at  the  bottom  of  the  ocean.  Curves  presented  there 

(2) 

compared  the  model  with  some  actual  measurements v  '  and  showed  good 
agreement  at  low  sea  state.  It  was  3hown  that  the  largest  amount  of  noise 
power  was  received  at  nearly  horizontal  angles,  and  that  the  noise 
intensity  decreased  as  the  vertical  angle  increased,  reaching  a  value 
10-20  db  lower  when  the  array  was  beamed  straight  up. 

The  question  arcse  as  to  whether  this  same  general  vertical 
distribution  of  noise  would  be  observed  by  a  hydrophone  array  near  the 
surface,  say  between  100  and  500  feet  below  the  surface.  This  paper 
undertakes  to  answer  that  question.  The  model  of  surface  noise  discussed 
by  Talham  is  used  as  is  the  general  fabric  of  his  derivation.  However, 
since  the  hydrophone  is  near  the  surface,  the  fine  structure  of  the 
sound  velocity  profile  near  the  surface  becomes  important,  and  consequently 
different  profiles  are  used.  The  one  employed  by  Talham  ignored  this  fine 
structure  near  the  surface  as  being  insignificant  to  measurements  made  at 
the  bottom  of  the  ocean. 
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I I .  Derivation  of  Vertical  Distribution 


The  geometry  used  by  Talham  with  slight  modifications  is  shown  in 
Figure  1.  A  hydrophone  located  at  depth  H  is  directed  at  an  angle  8  from 
the  horizontal,  and  for  simplicity  is  considered  to  be  omnidirectional  with 
respect  to  azimuth.  Noise  sources  are  assumed  to  be  located  uniformly 
over  the  surface  with  density  D  sources  per  unit  area.  The  sources  are 
statistically  independent,  and  each  emits  P»g(0^)  watts  per  unit  solid 
angle  into  the  water  at  an  angle  8^  with  respect  to  the  horizontal. 


Figure  1  Geometry  of  the  Model 

Talham  shows  that  the  intensity  (watts  per  unit  area)  of  noise 
arriving  at  the  hydrophone  between  the  angles  and  0q  +  A0  ,  and  due 
to  those  sources  lying  in  the  innermost  surface  ring  (ring  no.  l)  of 
area  As,  is 

AI1  »  DPg(81)W  exp  |-2aRjj  jsin  oj  ^2n  cos  OqA0J  (l) 

where  P  and  D  are  defined  above  and  W  is  the  ratio  of  backward-to-forward 
spreading  loss,  which  plays  a  significant  role  in  Talham 's  derivation  and 
is  given  by 
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W  =  W(0Q) 


(cos 
cos  e 

o 


(2) 


The  attenuation  of  the  noise  energy  is  given  by  the  exponential  function 
where  a  is  the  attenuation  constant  and  R^  is  the  ray  path  length  from 
the  hydrophone  to  the  surface.  The  factor  2n  cos  0qA6  is  the  solid 
angle  at  the  hydrophone  inside  of  which  noise  from  As  can  reach  the 
hydrophone.  The  factor  2n  cos  0q  accounts  for  the  fact  that  this  solid 
angle  decreases  as  Qq  becomes  larger.  The  inverse  sin  dependence 
takes  into  account  the  rapid  increase  in  As  as  6^  decreases,  and  therefore 
the  rapid  increase  in  the  number  of  noise  sources  contributing  to  AI. 

This  dependence  is  most  clearly  seen  in  the  following  argument,  where  for 
simplicity  the  rays  are  taken  as  straight  lines.  Thus  ■  0  ,  and 

as  shown  in  two  dimensions  in  Fig.  2,  the  areas  AA  and  As  are  related  by 
AA  «■  As  sin  for  small  AS. 


Figure  2  Simplified  Geometry  Illustrating  [sin  0.1  Dependence 


Ignoring  attenuation  and  considering  the  usual  inverse  square  law 

for  spreading  loss,  each  noise  source  in  As  produces  at  H  an  intensity 

2  2 
proportional  to  l/R^  .  However,  AA  is  proportional  to  cos  ©o  for 

2 

a  fixed  AO  and  consequently  As  is  proportional  to  R^  cos  0 ysin  . 


C-3 


2 

Thus  the  number  of  noise  sources  in  As  is  proportional  to  R^  cos  9 ^/sin  0^ 
and  the  aggregate  intensity  at  H  varies  as  cos  9 ^sin  9^  . 

Finally,  0^  and  R^  are  related  to  9q  by  standard  ray  tracing  relations. 
For  instance,  by  Snell's  law  we  have 


Chyd 
cos  9 

o 


Csurf 
cos  9^ 


(3) 


where  c^  ^  and  csu^  are  the  velocities  of  sound  at  the  hydrophone  and 
surface  respectively.  The  constant  is  peculiar  to  the  ray  arriving 
at  the  hydrophone  at  angle  9q  and  is  normally  called  the  'ray  parameter.' 

Also  contributing  to  the  total  noise  arriving  at  angle  Qq  is  energy 
that  arrives  at  As  along  rays  originating  at  greater  distances  on  the 
surface,  as  shown  in  Fig.  3  in  two  dimensions. 
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The  cases  of  angles  Qq  below  horizontal  are  also  shown  to  Indicate 

the  similarities  and  differences  between  the  two  cases.  Since  a  given 

ray  has  associated  with  it  a  particular  parameter  <y  by  (3 ),  at  each 

depth  all  the  rays  that  finally  reach  the  hydrophone  at  angle  tOQ 

exhibit  this  parameter,  and  thus  in  Fig.  3  all  rays  shown  have  angle  0^ 

as  illustrated,  where  c(y  ) 

cos  02  -  AO 

and  c(yg)  is  the  speed  of  sound  at  depth  y2>  (Note  that  since  y2  is 
positive  downward,  angles  are  measured  positive  in  a  clockwise  direction. ) 
This  symmetry,  which  requires  that  the  surface  and  bottom  be  parallel, 
simplifies  the  addition  of  ray  components.  We  now  consider  the  two  cases 
separately. 

III.  Upward-Looking  Case  (O0  <  0) 

From  Fig.  3  it  is  seen  that  there  exist  components  of  the  total 
intensity  at  angle  0Q  originating  at  x^,  Xy  etc.  The  ray  originating  at 
x^  travels  a  total  path  length  of  2R  +  «  3R^  +  23^  before  reaching  the 

hydrophone,  and  makes  a  reflection  off  both  the  bottom  (reflection 
coefficient  p )  and  the  surface  (reflection  coefficient  V ).  Since  the 
geometry  is  otherwise  the  same  as  for  AI^,  and  since  phase  shift  is 
unimportant  (the  noise  sources  are  independent),  the  contribution  due  to 
the  sources  at  is,  from  (l): 

AI3  ■  AljPr  exp[-luxR]  (5) 

Consequently  the  total  intensity  at  H  due  to  all  such  contributions  is 


l-  yp 


(6) 


To  obtain 
solid  angle  at 
through  by  AO, 
vertical  noise 


the  final  form,  we  remove  the  dependence  on  the  size  of  the 

the  hydrophone  by  deleting  2n  cos  0q  in  (1),  divide 

and  let  AO— *0  .  This  yields  for  the  (upward-looking) 

distribution  N  (0  ) 
up  o' 

D  Pg(0  )  W  exp(-2a  R,  ) 

N  (0  )  -  - - - - - 

up'  o  r 

(sin  1-  py  exp(-ljaR) 


Another  way  of  considering  the  deletion  of  2n  cos  is  to  consider  the 
hydrophone  array  to  be  no  longer  aziimthally  omnidirectional,  but  to  have 
an  azimuthal  beam  width  of,  3ay,  o,  measured  normal  to  the  ray  entering 
at  angle  0q.  This  replaces  2n  cos  Qq  in  (l)  by  the  single  term  a, 
which  then  of  course  is  retained  in  (7 ).  As  the  final  results  are  to  be 
normalized  anyway,  the  two  approaches  are  equivalent. 

It  is  possible  that  the  hydrophone  can  be  located  at  a  depth  such 
that  rays  arriving  at  angle  0Q  cannot  reach  the  bottom.  This  is  true 
whenever 

<  cbottom 


where  c^0^0m  is  the  speed  of  sound  at  the  ocean  bottom.  A  typical 
example  is  shown  in  Fig.  U  along  with  a  velocity  profile.  For  a  given  6q, 
Gy  given  by  (3)  is  a  sound  speed  which  is  equal  to  the  speed  of  sound  at, 
say,  depth  d„  Thus  by  (i;)  the  ray  angle  must  be  zero  at  depth  d.  Since 
cos  Og  in  (U)  cannot  exceed  unity,  no  ray  with  parameter  Gy  as  shown  can 
penetrate  beneath  depth  d.  Thus  in  ^Up(®0 )  of  (7 ),  the  term  p  is  removed 
(or  set  to  unity),  and  path  length  R  is  as  shown. 
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IV.  Downward-Looking  Case  ( 0  >0) 

Following  the  derivation  for  the  upward-looking  case,  it  is  clear  that 
the  only  difference  in  the  two  cases  is  that  the  ray  reaching  the  hydrophone 
at  an  angle  below  horizontal  must  have  travelled  further,  and  if  the 
bottom  is  reachable,  will  have  sustained  an  additional  bottom  reflection. 

Thus  from  Fig,  3  it  is  seen  that  a  ray  must  travel  further  by  a 
distance  R  +  Rg  -  R^  =  2R^  than  a  ray  arriving  above  the  horizontal, 
and  that 

Ndo,m<V- <?) 


Again,  if  the  bottom  cannot  be  reached,  p  is  set  equal  to  unity. 

V.  Various  Propagation  Situations 

As  was  suggested  by  the  sample  velocity  profile  given  in  Fig.  U  above, 
°V  "  chyc/cos  °o  deteri7d-nes  a  velocity  of  sound  for  each  angle  Qq.  Given 
a  particular  profile,  G^.  thus  delineates  various  regions  in  the  ocean  in 
which  the  given  ray  may  exist.  It  is  interesting  to  consider  the 
possibilities  as  a  function  of  0q. 
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A  typical  contour  is  given  in  Figure  $  along  with  two  different 
hydrophone  depths,  and  Hg.  Associated  with  are  three  rays  1,2,3 
showing  the  three  possibilities  for  rays  arriving  at  H^.  Ray  1  is 
confined  to  the  region  near  the  surface  because  <  c^^  ,  where 

cmode  k*18  speed  of  sound  at  the  peak  in  the  profile.  Ray  2  has  a  very 

long  path  length  for  all  contributions  but  the  nearest,  and  also  ray  2 

does  not  reach  the  bottom  since  3  is  seen  to  reach 

the  bottom,  as  will  all  rays  arriving  at  angles  greater  than  that  of  ray  3. 

Clearly  for  all  rays  do  reach  the  surface  and  thus  there  is  no 

null  in  the  intensity  due  to  surface  noise.  This  is  not  true  in  the  case 

of  a  hydrophone  at  Hg,  as  shown  by  ray  A,  which  cannot  reach  the  surface 

(or  originate  at  the  surface)  since  <  Cmode  .  Thus  a  null  in 

A 

surface  noise  is  expected  for  all  angles  0Q  satisfying 


cos  6  > 

o 


(null  exists) 


(10) 


mode 


where  c„  is  the  speed  of  sound  at  H0.  As  0  increases  further,  so  that 
2  ^  o 

this  condition  is  no  longer  satisfied,  then  rays  similar  to  ray  B  or 
ray  C  are  observed. 

Combining  these  thoughts  into  a  single  representation,  one  expects 

that  a  polar  plot  of  N  (0  )  and  N.  ,  (0  )  vs.  0  will  exhibit  fairly 

up  o  clown  o  o 

well-defined  regions,  each  corresponding  to  one  of  the  ray  types  discussed. 
This  is  suggested  in  Figure  6,  where  the  contour  and  depths  from  Figure  5 
have  been  used. 
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3.  deep  water  - 

2. 


1;  near  surface 

"1. 


3. 


hydrophone  at 

Figure  6  Polar  Delineation  of  Ray  Type  Regions 

It  is  expected  that  there  will  be  a  rather  rapid  decrease  in  N  ( 0q  ) 
and  N(jown(0o)  at  the  boundary  between  regions  1  and  2,  since  the  path 
length  abruptly  becomes  much  larger,  and  consequently  contributions  Aly 
Alt-,  etc.,  are  more  attenuated.  Such  a  rapid  decrease  is  not  expected  at 
the  boundary  separating  regions  2  and  3,  however,  since  of  necessity  the 
grazing  angle  with  the  bottom  is  very  small  near  this  boundary,  and 
consequently  the  reflection  coefficient  is  near  unity. 

For  completeness  it  should  be  remarked  that  for  profiles  where 

cbottom  <  Cmode  *  re8ions  2  and  0  are  absent,  since  then  any  ray 
penetrating  into  deep  water  will  always  reach  the  bottom. 

VI.  Computation  of  Noise  Distribution 

It  is  noted  from  (7)  that  the  vertical  distribution  \p(®0)  [■and 

consequently  Ndown(0o}]  re<luires  that  °r  h>  R>  and  P  =  P(0bottom) 
be  known.  To  this  end  a  ray  tracing  program  was  written,  and  this  was 

applied  to  several  profiles  approximated  piece-wise  by  straight  lines. 
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The  ray  tracing  program  is  described  in  the  Appendix.  The  significant 
outputs  of  this  program  are  and  Rg,  the  ray  path  lengths  from  the 
hydrophone  to  the  surface  and  bottom,  respectively,  p  was  found  using 

(3) 

the  curves  shown  in  Figure  7,  originally  reported  by  Marsh  and  Schulkin. 

The  same  surface  reflection  coefficient  y  and  attenuation  constant  a  as 

used  by  Talham  were  also  used,  namely,  y  -  1  ,  and  c  -  .0038  f^^  Np/kyd  , 

where  f  is  in  kilocycles  per  second.  A  different  value  of  y  is  also  tried, 

as  discussed  below.  Two  source  directionality  functions  g(01)  were  used, 

g^)  =  1  and  gCft^)  -  sin  01  ,  although  others  were  checked  and  showed 

no  peculiar  features.  Talham  showed  that  at  least  on  the  ocean  bottom, 

and  at  low  sea  states,  g(Q^)  »  1  agreed  most  closely  with  available  data. 

The  results  obtained  are  presented  in  Figures  10-18  as  polar  plots 

of  N  (0  )  and  N , _ (0  ).  All  plots  are  in  db  relative  to  unity,  with 

up  o  clown  o 

P  and  D  normalized  to  one.  Two  values  of  frequency  were  used,  U00  cps 
and  1000  cps,  as  being  representative  of  the  range  of  interest.  Four 
sound  velocity  profiles  ^  were  used,  as  given  in  Fig.  8.  Attention  was 
concentrated  on  the  Bermuda  profiles,  while  the  Iceland  profiles  are 
included  for  purposes  of  comparison.  Various  features  of  these  graphs 
are  discussed  below. 

Certain  situations  can  be  handled  analytically  for  comparison  with 
computer  results,  and  to  show  the  major  effects  contributing  to  NUp(®0)* 

In  particular  the  cases  where  0Q  is  near  0°  or  90°  can  be  studied. 

Assume  that  ®0  “  f  •  Then  from  (7 )  we  have,  since  =  ®0  “  ■§  » 
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(11) 


DPg(f)e-2aH 

l-Vp(f)eJ,a5 


DPg(f) 


(12) 


where  H  and  B  are  the  depths  of  the  hydrophone  and  bottom  respectively,  and 
the  dependence  of  0  on  the  bottom  angle  is  indicated  as  p(^)  .  The 
approximation  in  (12 )  stems  from  the  fact  that  at  frequencies  in  the  hundreds 
of  cycles  the  exponents  in  (11)  are  very  small. 

At  the  other  extreme,  when  the  hydrophone  lies  in  a  region  of  positive 
gradient,  as  at  in  Fig.  5,  and  0Q  is  small  enough  so  that  the  ray  is 
confined  to  depths  less  than  ^  (region  1  of  Fig,  6),  then  and  Rg 
are  reasonably  small,  W  1  ,  and  (5  -  1  .  Thus  from  (7)  NUp(^0 )  becomes 
approximately,  with  Y  -  1  , 


V°o> 


Dpg(a, ) 


sin  Q^(IkxR) 


(©o  small) 


(13) 


This  may  be  reduced  further  if  the  velocity  gradient  0  at  the 
hydrophone  is  approximately  constant  up  to  the  surface.  Then  the  radius 
of  curvature  r  of  the  ray  may  be  shown  ^  ^  to  be  r  ■  G^/G  .  Using  the 
geometry  of  Fig.  9,  where  Rg  is  the  distance  from  H  to  the  maximum  depth 
reached,  it  is  seen  that 


and 


\  ■  r<®1-00) 


R3  ‘  r0o 


cos  G 

0 


COS  G^ 


H 

r 


(1U) 

(15) 
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But  since  C y  -  0hy</C0S  ®0  *  (15)  yield8 

cos  ©^  -  (1  —  k )  cos  0^ 

where 

,  ,  ,  HG  csurf 

chyd  chyd 


(16) 

(17) 


Thus  (13 )  now  becomes 


DPg(e, ) 

N  (o  )  -  - i - 

UP  0  i*a  sin  ©^  r  ©^ 


(18) 


Using  cos  0  -  1  -  -y  twice  in  (16),  (18)  then  yields 


N  (0  )  = 
up'  o' 


DPg(©_  )  G  cos  0 
1  o 


Ua  sin  01  Chy(J  ^2k  *  (1-  k)  Qq2 


(0  small ) 
0 


The  two  forms  for  g(01)  treated  here  are  g^)  *  1  and 
g(Q^)  ■  sin  ©^  .  These  forms  yield  the  results 


DP  G  cos  © 

N  (0  )  -  0 

up'  o 


*‘av[2kt(1-lc)  °o2. 


if  g(e1)  *  l 


(19) 


(20) 


D  PG  cos  © 

Veo> - -  if  g(e1).Si„a1 

chyiTk*  (1"k)  V  (2. 

The  approximation  sin  ©1  =  is  used  in  (20).  As  an  example,  for 

the  Bermuda  winter  profile,  with  the  hydrophone  located  at  100  ft,  k  has 

o 

the  value  .0003  .  Thus  the  term  2k  is  less  than  (l-k)O  if  0  >  1.4 

o  o  — 

degrees.  The  forms  of  (20)  and  (21)  show  that  as  ©Q — >0  ,  Nup(©0) 
approaches  the  values 
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(0) 


DP 

Hal 

TOfVsr 


if  g(ex)  =  i 
if  g(01)  =  sin  0^^ 


(22) 


It  is  noted  that  the  only  dependence  on  the  velocity  profile  for 
0Q  *  0  occurs  in  k.  Thus,  for  example,  in  any  ocean  for  which  a  constant 
positive  velocity  gradient  exists  from  the  surface  to  the  hydrophone,  and 
for  H  =  100  ft  and  f  *  U00  cps  ,  we  have  in  db  (with  P  =  D  =  1  ) 


35.9  db 


V0)! 


db 


35.9  -  5  log10  2§r  db 


if  g(e1)  -  l 
if  g(01)  =  sin  01 


(23) 


For  both  the  Bermuda  and  Iceland  winter  profiles  (23)  yields  19. 6  db  for 

the  g(O^)  =  sin  6^  case.  These  agree  closely  with  the  computer  results, 

as  seen  in  Figures  10  and  15.  When  the  hydrophone  lies  in  an  ocean  region 

for  which  the  gradient  is  negative,  as  at  Hg  in  Fig.  5,  the  simple 

geometric  techniques  used  above  fail  due  to  the  complicated  ray  path. 

Shadow  regions  are  easily  found,  as  mentioned  in  Section  17,  but  actual 

results  for  N  (0  )  require  a  computer.  In  all  cases  N.  (0  )  is  easily 
up  o  down  o 

found  from  N  (0  )  by  means  of  (9). 

UP  ® 

It  is  instructive  to  consider  the  physical  causes  for  the  various 
results  obtained  in  (22 ).  Beginning  with  (13 ),  it  is  seen  that  there  are 
two  separate  contributions  to  ^Up(®0  )•  The  first  is  the  inverse  sin  0^ 
effect,  which  arises  geometrically  (see  Fig.  2)  and  relates  the  number  of 
noise  source  elements  spanned  by  the  beam  of  the  hydrophone.  This  effect 
is  cancelled,  of  course,  if  it  is  assumed  that  g(Q^)  -  sin  0^.  Thus  the 
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magnitude  of  this  effect  can  be  judged  from  a  comparison  of  the  graphs  with 

g(01)  *  1  and  g(©1)  -  sin  ,  as  in  Figs.  10,  11,  15,  and  16.  The 

effect  is  present  in  the  other  cases  as  well,  of  course,  but  (22)  is  only 

applicable  to  the  cases  corresponding  to  these  figures.  The  second  effect 

is  that  of  the  inverse  path  length  dependence,  which  is  easily  traced  to  (6) 

and  is  due  to  the  summation  of  contributions  from  remote  surface  elements. 

It  is  thus  seen  that  these  remote  sources  add  up  significantly  at  small 

values  of  Qq,  producing  as  in  (18)  an  inverse  0^  dependence  (where  0^  is 

small  when  Qq  is  small).  This  is  true  for  low  frequencies  since 

attenuation  is  then  small.  The  magnitude  of  this  effect  is  demonstrated 

in  Fig.  lU,  where  the  complete  N  (Q  )  and  N .  (0  )  are  compared  with  the 

up  o  down  o 

corresponding  ambient  noise  levels  when  only  the  first  contribution  to  the 
noise,  that  from  ring  no.  1  (see  Fig.  1),  is  considered.  These  values 
were  obtained  by  setting  V  =  0  in  (6).  It  is  seen  that  the  noise  level 
due  to  the  first  contribution  is  independent  of  frequency  for  upward- 
looking  cases  and  only  slightly  dependent  for  downward-looking  cases. 

The  independence  is  due  to  the  negligible  attenuation  for  the  short  path 
lengths  involved,  while  the  slight  dependence  in  the  downward  direction  is 
due  to  the  fact  that  the  bottom  reflection  coefficient  p  is  a  decreasing 
function  of  frequency  [Fig,  7J. 

This  second  effect  due  to  contributions  from  remote  sources  is 
heavily  dependent  on  the  fact  that  the  surface  reflection  coefficient  V 
is  assumed  to  be  equal  to  unity.  If  this  is  not  the  case,  but  V  differs 
from  unity  by  6,  (y  ■  1-  6),  then  (7)  leads  to 
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DPg(0,  ) 

N  (Q  )  *  - s - = - r 

up  0  sin  01[l-y(l-UaR)l 

DPg(Q1) 
sin  9^  [6+  liaHy] 


(O0  small)  (22*) 

(25) 


Marsh  et  al.  have  derived  a  theoretical  dependence  for  y  on  several 
parameters,  ^  arriving  at 

6  *  .1*85  f1*5  (1.77  h)1'6  sin  91  (26) 


where  f  is  the  frequency  in  kcps  and  h  is  the  rms  wave  height.  This 
dependence  is  valid  for  small  9^  and  low  sea  state.  However,  it  will  have 
little  effect  on  NUp(®0)»  since  0^  is  very  small  when  0o  is  near  zero,  and 
consequently  y  «  1  .  A  sample  computer  result  is  shown  in  Fig.  10,  which 
compares  the  values  of  Nup(°0)  when  y  «*  1  and  when  y  **  1-  6  .  It  is  seen 
to  cause  a  negligible  reduction  in  ambient  noise  level. 

VII.  Further  Discussion  of  Results 

The  rather  peculiar  shape  of  the  vertical  noise  distributions  of 
Figs.  10-18  merits  some  discussion.  Examining  Fig.  10  as  an  example,  it 
is  seen  that  the  noise  level  drops  suddenly  as  Qq  passes  beyond  U°  or  5°. 
It  is  seen  that  thi3  is  true  for  both  the  g(O^)  ■  1  and  g(O^)  ■  sin  0^ 
cases,  so  that  the  effect  cannot  be  due  to  the  falling  off  of  the  number 
of  noise  elements  involved  (the  inverse  sin  0^  effect).  The  effect 
instead  is  that  of  attenuation,  for  as  shown  in  the  table  of  R^  and  Rg 
accompanying  Fig.  10  (which  applies  as  well  to  Fig.  11),  there  is  a 
sudden  jump  in  R^  +  Rg  as  ©q  reaches  5°,  for  a  hydrophone  at  100  ft. 

This  is  due  to  a  jump  of  the  ray  from  a  region  1  type  to  a  region  2  type, 
as  shown  in  Fig.  6.  The  angle  0Q  ■  5°  is  just  sufficient  to  allow  the 
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ray  to  pass  out  of  the  shallow  water  channel  into  a  deep  water  path.  This 
causes  all  contributions  but  the  first  to  the  total  noise  level  to  travel 
much  greater  distances  to  the  hydrophone.  Due  to  the  attenuation  of  sound 
in  water,  these  now  very  remote  contributions  are  much  smaller,  and  the 
total  noise  level  is  consequently  reduced. 

Beyond  this  sudden  drop  there  is  a  gradual  decay  in  the  noise, 
except  in  the  case  where  g(O^)  -  sin  6^  ,  where  slight  increases  are 
seen  due  to  the  fact  that  +  Rg  again  decreases  slightly.  The  decay 
increases  as  soon  as  the  bottom  is  reached  (which  occurs  at  12°  in  all 
curves  of  Figs.  10  and  11),  for  the  bottom  reflection  coefficient  begins 
to  decrease  from  unity.  Finally  this  levels  out  as  seen  from  P-ig.  7  and 
thus  the  g(Q^ )  -  sin  ©1  case  is  essentially  constant  while  the  g(®1)  3  1 
case  decreases  due  to  the  inverse  sin  effect. 

The  difference  in  hydrophone  depths  is  seen  to  be  unimportant  except 

at  small  values  of  since  at  small  SQ  the  slightly  increased  path 

length  for  larger  hydrophone  depths  permits  a  small  change  in  0^  and 

N  (0  )  is  very  sensitive  to  0.  at  small  values  of  0  .  At  larger  values 
up  O  X  o 

of  0Q  the  rays  are  geometrically  very  similar  for  different  hydrophone 
depths. 

In  the  summer  profile  cases  of  Figs.  12  and  13  the  g(fi^)  “  sin  0^ 
case  shows  a  constant  noise  level  over  a  range  of  6q  near  n/2 .  This  is 
due  to  the  fact  that  R^  and  R^  +  Rg  are  essentially  constant  over  this 
range,  as  a  sketch  of  ray  paths  for  such  a  profile  will  readily  show. 

Since  the  inverse  sin  6^  effect  is  cancelled  in  this  case,  the  constancy  of 
the  noise  is  understandable.  It  ends  as  soon  as  the  bottom  is  encountered 
(at  8°  for  a  hydrophone  depth  of  100  ft,  and  at  11°  for  a  depth  of  500  ft). 
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For  still  larger  angles  the  features  mentioned  above  again  apply.  The 
summer  cases  also  exhibit  regions  in  which  a  null  exists  in  the  ambient 
noise.  These  are  considered  in  more  detail  below. 

VIII.  Conclusions 

It  has  been  shown  that  an  ambient  noise  model  which  assumes  a 

uniform  surface  distribution  of  independent  noise  source  elements  gives 

rise  to  one  of  two  very  different  noise  distributions  vs.  vertical  angle. 

If  the  hydrophone  lies  at  a  depth  where  a  negative  velocity  gradient 

exists,  there  will  be  no  surface  ambient  noise  received  at  angles  near 

the  horizontal,  although  at  slightly  larger  angles  the  noise  intensity 

will  rapidly  jump  to  a  high  level.  If  the  hydrophone  lies  in  a  region 

of  positive  velocity  gradient  (and  if  the  sound  velocity  at  the  hydrophone 

is  greater  than  that  at  the  surface),  then  no  shadow  region  exists  and 

the  noise  level  grows  rapidly  as  the  vertical  angle  approaches  the 

horizontal.  This  peaking  near  0Q  *  0  is  due  to  two  important  features. 

The  first  occurs  because  at  small  values  of  0  a  small  solid  angle  at 

o 

the  hydrophone  covers  a  very  large  area  of  surface  due  only  to  the 
geometry  involved.  The  second  arises  from  contributions  to  the  noise 
generated  far  from  the  hydrophone  and  propagating  to  the  hydrophone  via 
many  surface  bounces.  This  effect  may  be  greatly  reduced  if  the  surface 
reflection  coefficient  is  heavily  dependent  on  sea  state.  The  first 
effect  may  be  counteracted  partially  if  the  surface  noise  elements  are 
directional,  such  that  more  power  is  transmitted  at  steep  downward  angles 
than  at  very  shallow  angles  to  the  surface. 

The  results  in  the  winter  profile  cases  shown  above  are  seen  to 
depend  only  to  a  small  extent  on  frequency  and  hydrophone  depth,  for  low 
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frequencies  (100-1000  cps )  and  shallow  depths  (100-500  ft).  Except  for 
the  seasonal  shifts  in  velocity  profiles  which  cause  radical  changes  in 
the  noise  distribution,  the  distribution  is  rather  insensitive  to  profile, 
as  seen  by  comparing  the  Bermuda  and  Iceland  cases  in  each  season. 

The  superficially  promising  null  in  the  noise  at  small  values  of  Qq 
during  the  summer  months  is  not  actually  very  helpful,  since  it  is  due  to 
sound  channeling  effects.  Thus  a  target  near  the  surface  would  also  be 
hidden  from  the  hydrophone.  The  exception  to  this  would  occur  if  the 
submarine  and  hydrophone  were  both  in  the  negative  gradient  region  (below 
Dm  ode  in  Fig.  5)  and  at  approximately  the  same  depth.  Then  the  submarine 
signal  would  arrive  at  the  hydrophone  at  an  angle  near  zero,  for  which 
there  is  a  null  in  the  surface  ambient  noise.  This  case  possibly  could 
be  taken  advantage  of,  in  order  to  combat  such  noise. 

The  results  indicate  that  bottom  bounce  techniques  might  be  employed 
in  order  to  take  advantage  of  the  greatly  reduced  ambient  noise  at  angles 
steeper  than  about  20°.  The  noise  level  at  these  angles  can  be  as  much 
as  28  db  below  that  at  near  horizontal  angles. 

The  above  results  depend  strongly  on  the  assumed  noise  model.  The 
need  for  careful  measurements  of  the  vertical  noise  distribution  near 
the  surface  therefore  appears  to  be  indicated.  By  comparing  such  results 
with  the  various  cases  given  in  Figs.  10-18,  it  should  be  possible  to 
make  some  judgments  as  to  the  functional  form  of  g(Q^),  the  surface 
reflection  coefficient  y,  and  indeed  the  validity  of  the  basic  model 
concept,  that  of  a  uniform  surface  distribution  of  noise  sources. 
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Appendix  The  Ray  Tracing  Program 

In  order  to  evaluate  NUp(®0)  and  Ndown^9o^  following  quantities 
must  be  computed. 

0^s  the  angle  of  intersection  of  the  ray  with  the  surface 
©bottom5  an€^-e  intersection  of  the  ray  with  the  bottom 
R^i  the  ray  path  length  from  the  hydrophone  to  the  surface 
Rg*  the  ray  path  length  from  the  hydrophone  to  the  bottom 


If  the  ray  does  not  reach  the  bottom,  then  is  the  path  length  of 
the  ray  from  the  hydrophone  to  its  position  of  maximum  depth,  at  which 
point,  of  course,  it  is  moving  horizontally. 

If  the  velocity  profile  is  approximated  by  several  straight  line 
segments,  the  ray  path  in  the  region  of  depth  for  each  segment  is  uniquely 
determined  by  the  angle  with  which  it  enters  the  region,  the  velocity  of 
sound  at  the  depth  at  which  it  enters,  and  the  velocity  gradient  in  the 
region. 

By  Eq.  (h)  the  angle  at  each  depth  is  known  as  soon  as  ©q  and 
are  given,  and  the  sound  velocity  at  that  depth  is  computed.  Thus  ©^ 
and  ®bottom  (when  applicable)  can  be  found  without  any  ray  tracing  as  such. 

The  program  used  begins  at  the  surface,  and  taking  each  linear 
segment,  subdivides  the  corresponding  depth  interval  into  a  sufficient 
number  of  subintervals.  The  criterion  used  was  that  the  ray  angle  should 
not  change  by  more  than  a  specified  fraction  of  a  degree  over  the  subinterval. 
A  change  of  one  degree  was  found  tolerable.  Thus  in  this  subinterval  the 
ray  path  is  essentially  straight,  and  if  it  is  directed  at  an  angle  ©^ 
with  respect  to  the  horizontal,  the  path  length  in  the  subinterval  is 


032 


d^/sin  9^  ,  -where  d^  is  the  change  in  depth  in  the  subinterval.  Thus 
is  merely  the  sum  of  all  these  contributions,  until  the  depth  of  the 
hydrophone  is  reached.  This  number  is  stored  and  then  the  computation 
resumes  downward  to  find  R^.  When  either  the  bottom  or  an  angle  -  0 
is  reached,  the  computation  is  terminated. 
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Summary 

The  detection  of  s  single  sonar  pulse  by  means  of  a  correlation  detector 
is  anal  zed  in  the  absence  of  Doppler  shift  due  to  transmitter,  target,  or 
scntterer  motion.  The  returned  signal  is  assumed  to  be  a  delayed  replica  of 
the  transmitted  signal.  Noise  fields  consisting  primarily  of  volume  rever¬ 
beration,  surface  reverberation,  and  embient  noise  are  considered  separately. 

The  reverberation  noise  is  assumed  to  be  generated  by  independently  located, 
Poisson  distributed  scatterers,  dispersed  throughout  the  illuminated  volume 
(volume  reverberation)  or  near  the  illuminated  surface  (surface  reverberation). 

The  transmitted  signal  is  assumed  to  be  a  pulse  of  sinusoidal  carrier  with 
or  without  superimposed  linear  frequency  modulation.  The  following  results 
are  obtained: 

1)  The  bandwidth  of  the  transmitted  signal  imposes  an  effective  range 
gate  on  output  noise  arising  from  the  intermodulation  products  of 
signal  with  reverberation.  The  width  of  this  gate  is  of  the  order 
of  the  velocity  of  sound  divided  by  the  bandwidth  of  the  signal. 

(  this  ratio  may  be  regarded  as  the  corrclr.tion  distance  in  the 
water  of  the  transmitted  signal).  Only  scatterers  whose  range  differs 
from  that  of  the  target  by  no  more  than  the  gate  width  contribute  to 
this  type  of  output  noise  (which  is  dominant  under  high  signal  to 
noise  conditions). 

2)  Output  noise  consisting  of  the  intermodulation  products  of  reverberation 
with  reverberation  is  contributed  by  all  scatterers  illuminated 
simulaiiuously  with  the  target.  Contributing  scatterers  therefore 

lie  within  a  range  gat^.  determined  by  the  pulse  duration  rather 
than  the  (possibly  much  shorter)  signal  correlation  time.  However, 
the  two  members  of  each  contributing  pair  must  be  separated  form  each 


other  by  a  range  difference  no  greater  than  the  "range  gate"  defined 
in  1).  Intermodulation  of  reverberation  with  reveration  is  the  dominant 
noise  at  low  input  signal  to  reverberation  ratios.  Except  for  the 
phenomenon  discussed  in  3),  there  is  no  angular  discrimination 
against  either  type  of  reverberation  noise  aside  from  the  obvious 
effects  of  transmitter  and  receiver  beam  patterns. 

3)  A  special  case  of  2)  is  output  noise  due  to  the  intermodulation  of 
the  returns  from  a  particular  scatterer  to  the  two  receivers.  This 
type  of  noise  is  important  only  when  one  or  a  few  strong  scatterers 
(false  targets)  contribute  a  significant  fraction  of  the  total 
reverberation.  False  target  returns  are  subject  to  angular  dis¬ 
crimination  .being  sharply  attenuated  if  they  are  separated  from 
the  target  by  more  than  a  certain  angle.  This  angle  defines  the 
minimum  angular  resolution  of  which  the  detector  is  capable.  For 
a  broadside  target  it  is  the  arc  sine  of  the  ratio  of  the  signal 
correlation  distance  to  the  spacing  between  receivers.  The  resolution 
is  best  in  the  broadside  direction  and  poorest  in  the  endfire 
direction. 

li)  In  the  absence  of  a  target  the  average  detector  output  is  zero  except 
for  a  transmitter-receiver  combination  with  very  narrow  beam  patterns 
trained  very  close  to  the  axis  of  the  receiving  amy.In  the  latter 
case  a  non-zero  average  output  may  result  from  the  inability  of 
the  detector  to  resolve  different  scatterer  groups  within  the  illuminate 
area. 

5)  In  the  absence  of  strong  scatterers  (false  targets)  and  for  beams 
trained  far  enough  from  the  array  axis  so -that  the  resolution 
problem  mentioned  in  4)  does  not  arise,  the  reverberation  may  be 
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regarded  as  a  Gaussian  noise  process.  For  a  fixed  input  signal  to  noisi 
ratio,  there  is  then  no  difference  in  performance  between  a  detector 
operating  in  a  noise  field  consisting  primarily  of  reverberation  and 
a  similar  detector  operating  in  an  ambient  n6ise  field  with  the 
same  power  spectrum.  If  the  ambient  noise  field  is  broad-band 
but  is  filtered  so  that  its  spectrum  comes  close  to  matching  that 
of  the  signal,  there  is  still  no  important  difference  between 
detector  performance  in  reverberation  limited  and  ambient  noise 
limited  environments  (for  fixed  input  signal  to  noise  ratio).  One 
must  keep  in  mind,  however,  that  the  effective  input  signal  to  noise 
ratio  tends  to  vary  inversely  with  the  signal  bandwidth  for  ambient 
noise,  whereas  it  is  independent  of  the  signal  bandwidth  for  re¬ 
verberation.  Thus  maximum  output  signal  to  noise  ratio  is  generally 
achieved  by  the  use  of  broadband  signals  if  the  primary  problem  is 
reverberation  and  by  the  use  of  narrowband  signals  if  the  primary 
problem  is  ambient  noise.  In  the  former  case  one  takes  advantage 
of  the  range  gating  effect  produced  by  short  signal  correlation 
times,  in  the  latter  case  one  relies  on  noise  reduction  through 
the  use  of  narrow  filters  matched  to  the  signal  spectrum. 

6)  Under  the  conditions  stated  in  5)  the  output  signal  to  noise 

ratio  varies  as  the  square  of  the  input  signal  to  noise  ratio  for 
small  input  signal  to  noise  ratios  and  as  the  first  power  of  the 
input  signal  to  noise  ratio  for  large  input  signal  to  noise  ratios. 

Thus  the  correlation  detector  operates  much  like  a  coherent  detector 
for  large  input  signal  to  noise  ratios  but  more  like  an  incoherent 
detector  for  small  input  signal  to  noise  ratios.  The  output  signal 
to  noise  ratio  of  a  true  coherent  detector  (correlating  the  output 
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of  one  receiver  with  a  delayed  replica  of  the  transmitted  signal) 
is  precisely  twice  that  of  the  correlation  detector  in  the  limit  of 
high  input  signal  to  noise  ratio.  For  one  set  of  reasonable  parameter 
values,  the  correlation  detector  requires  an  input  signal  to  noise  ratio 
8  db  higher  than  that  of  the  coherent  detector  in  order  to  achieve  an 
output  signal  to  noise  ratio  of  6  db. 

7)  The  primary  difference  between  surface  and  volume  reverberation  is  the 
time  dependence  of  the  former,  caused  by  the  fact  that  surface  rever¬ 
beration  is  a  problem  only  when  elements  of  the  surface  ere  illuminated 
simultaneously  with  the  target.  This  leads  to  limitations  on  desirable 
pulse  length  which  are  not  present  in  the  case  of  volume  reverberation. 

8)  Reverberation  power  varies  with  range:  with  the  inverse  second  power  of 
range  for  volume  reverberation  and  the  inverse  third  power  of  range 

for  surface  reverberation.  The  input  signal  power  varies  with  the  inverse 
fourth  power  of  range.  Hence,  for  small  input  signal  to  noise  ratios, 
the  output  signal  to  noise  ratio  varies  as  the  inverse  fourth  power 
of  range  for  volume  reverberation  and  the  inverse  second  power  of  range 
for  surface  reverberation.  When  ambient  noise  dominates,  the  output 
signal  to  noise  ratio  varies  with  the  inverse  eighth  power  of  range. 

Under  conditions  of  high  input  signal  to  noise  ratio,  the  output  signal 
to  noise  ratio  varies  with  the  inverse  second  power  of  range  for  volume 
reverberation,  with  the  inverse  first  power  of  range  for  surface 
reverberation  and  with  the  inverse  fourth  power  of  range  for  ambient  noise. 
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I.  Detection  in  a  Reverberation  Limited  Environment 

This  report  is  concerned  with  the  detection  of  an  active  sonar  signal 
in  a  noise  environment  that  may  be  dominated  either  by  reverberation  or  by 
ambient  noise.  The  postulated  signal  processing  scheme  is  outlined  in  Figure  1. 
A  sonar  pulse  produced  by  a  transmitter  with  known  beam  pattern  is  reflected 
by  the  target  and  received  at  two  locations,  labelled  receiver  a  and  receiver  b 


Figure  1 

respectively.  Whether  each  of  these  receivers  consists  of  one  hydrophone 
or  of  a  group  of  hydrophones  is  immaterial  to  the  analysis  except  for  the 
effect  jr  the  signal  strength  received  from  any  given  direction.  As  long  as 
the  two  receivers  are  similar  and  the  target  is  relatively  remote  (rQ»>d),  any 
beam-forming  effects  of  the  receiver  can  clearly  be  lumped  with  those  of  the 
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transmitter,  riio  out.j>n+.  of  fine  receiver  is  delayed  by  an  amount  tf, 

r>K •*»on  so  as  to  bring  the  two  target  returns  into  alignment  (steering). 

It  is  then  multiplied  by  the  output  of  the  second  receiver  and  the  product 

y  is  smoothed  by  means  of  a  low  pass  filter. ^  Reverberation  is  assumed 

to  be  generated  by  a  series  of  scatterers  randomly  distributed  in  volume 

and  at  the  surface.  The  scatterers,  as  wall  as  the  source  and  target, 

2 

are  assumed  to  be  at  rest. 

The  tran5<mitted  signal  is  taken  as  a  linearly  frequency  modulated 
pulse  with  Gaussian  envelope.  Thus 


Transmitted  signal  =  exp 


cos(ea  t 
o 


(1) 


In  cases  of  practical  interest  the  pulse  width  Op,  maximum  frequency 
deviation  in  the  pulse  Kokand  mean  carrier  frequency  n satisfy  the 
inequalities 


\«  k»t  «  “o 


(2) 


In  a  reverberation  dominated  environment  the  output  v  of  receiver  a 


^Alternatively  or.e  might  use  transducers  a  and  b  for  transmission  as 
well  as  reception.  With  proper  delays  in  transmitted  and  received  signals 
the  signal  component  of  the  output  would  be  identical  with  that  obtained 
from  the  configuration  of  Fig.  1.  The  formal  expression  for  the  reverberation 
noise  would  be  more  cumbersome,  since  for  most  scatterers  the  round  trip 
times  to  one  receiver  from  the  two  sources  would  be  different.  It  is 
not  felt  that  the  inclusion  of  this  added  complexity  would  yield  significant 
additional  insight  into  the  detection  process, 

2 

For  the  relatively  wideband  signals  considered  in  much  of  this 
report,  the  Doppler  shifts  due  to  scatterer  motion  (even  due  to  target 
motion)  should  be  quite  small  compared  with  the  signal  bandwidth. 
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can  now  be  written  as  follows: 


V  (t)  =  — *•  exp 

V 


(t-ta) 

— r 


2- 

=os[%(t-ta)  +  £  (t-ta)2] 


'Zb 


exp  h 


1  i 


2i 


cos[coo(t-t1)  |  (t,-ti  )2  ] 


(3) 


r  +  r 
o  a 

In  Equation  (3)  t  =  — - -  is  the  round  trip  time  from  the 
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transmitter  via  the  target  to  receiver  a.^  c  is  the  velocity  of  sound 

in  water,  t^  is  the  round  trip  time  from  the  transmitter  to  scatterer 

i  and  back  to  receiver  a.  The  constant  A  depends  on  the  scattering 

cross-section  of  the  target  and  on  the  amount  of  radiation  received  by 

the  target,  i.e,  on  the  beam  parameters  of  the  transmitter.  In  similar 

fashion,  a^  depends  on  the  scattering  cross-section  of  the  i^  scatterer 

and  on  the  transmitter  beam  parameters.  The  dependence  of  the  signal 

and  noise  components  of  V  (t)  on  1/t  and  l/t.  respectively  reflects 

3  3  2. 

the  usual  spherical  spreading  loss  and  results  in  a  non-stationarity  of 
the  reverberation  noise. 

In  completely  analogous  fashion  the  output  7,  of  receiver  b  assumes 

D 

the  form; 


(h) 


Here  t^  and  are  the  signal  round  trip  times  from  the  transmitter  to 

_  j 

Strictly  speaking,  the  first  term  is  proportional  to  c  / (rora> .  With 
r  «r  this  becomes  approximately  4/t  ^ 
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receiver  h  via  the  target  and  scatterer  j  respectively.  V^(t)  is  delayed 
by  t  seconds  and  multiplied  by  V  (t)  to  yield  y(t).  After  some  algebraic 
simplification  one  obtains*. 
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If  the  weighting  function  of  the  low  pass  filter  is  w(t),  its 
output  z(t)  ic  given  by 


Sinusoidal  terms  in  2<u  t  have  been  discarded  since  thy  will  be 
filtered  out  by  the  low  pass°filter  in  any  case. 
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z(t)  = 


dp  w (p)  y(t-p) 


(6) 


For  computational  convenience  the  weighting  function  will  be  ehosen  as 

\2  n 


w(t) 


v7tc 


exp 


(7) 


This  corresponds  to  the  frequency  function 


G  (co)  =>  w(t)  exp  j-Jcfcj  dt 


exp 


aF  2 

~T  ^ 


expj-  joitD]  (8) 


As  long  as  the  time  delay  t  satisfies  t  »  ct^.  the  filter  is 

D  D  * 

approximately  realizable. 

Substituting  Equations  (5)  and  (7)  into  Equation  (6)  one  obtains 
after  extensive  algebraic  manipulations: 
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Equation  (9)  indicates  that  the  delay  ty  necessary  to  maximize 
the  signal  component  of  z(t)  is 


t  =  t  -t, 
r  a  b 


(11) 


With  this  adjustment  of  delay 
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{[(V‘a)-<Trtb)] [v K  -r  <*■  -  VYVtb  -  H']}  (12) 


The  signal  component  of  z(t)  clearly  peaks  at  t  =  t  +  tD.  The  statistical 

properties  of  the  random  variable  z(t  +  tD)  are  therefore  of  primary  importance 

in  any  discussion  of  detector  performance.^-  The  general  expression  for 

z(t  +  t  )  is  from  Equation  (12) 
a  d 


*(t+t  ) 


*  2  o 
A  c 


3  D'  2t  2t.  2  ffF 

3  u 


“^The  statement  that  the  output  is  observed  at  time  tfl+t  implies  selection 

of  the  proper  range  cell,  just  as  the  choice  t  =  t  -t.  implies  selection  of 

r  aD 

the  proper  bearing  cell.  If  bearing  and  range  are  not  known  a  priori,  one 
must  clearly  examine  all  possible  cells, 
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(13) 


For  ease  in  later  computations  it  will  be  convenient  to  distinguish 
between  two  cases: 

A )  The  bandwidth  of  the  transmitted  signal  is  determined  primarily 
by  the  frequency  modulation  (Ko^1»l/oT). 

B)  The  bandwidth  of  the  transmitted  signal  is  determined  by  the  pulse 
duration  (no  frequency  modulation,  K  =  0). 

In  case  A)  Equation  (13)  reduces  to  the  expression 


.  2  o 
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•  co3|[(ti.ta)-(Trtb)][vf^2  (yyt^)] 
The  equivalent  expression  for  case  B)  is 
a(t  +t.  ) 


(1U) 
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where  CT  -  - g - 2— 

°F  +  aT 


(15) 

(16) 


The  fird,  term  in  each  of  the  two  equations  ^(lU)  and  (l5)Jis  the  signal 
component,  the  next  two  terms  are  due  to  the  intermodulation  of  signal  from 
one  receiver  and  reverberation  from  the  other,  while  the  last  term  results 
from  the  intermodulation  of  the  reverberation  components. 

Several  qualitative  conclusions  concerning  detector  performance  can  be 
drawn  immediately  from  inspection  of  Equations  (lh)  and  (l5)s 

1)  If  the  signal  is  strong  compared  with  the  reverberation,  the  noise 
consists  primarily  of  intermodulation  products  of  signal  with 
reverberation,  and  the  last  noise  term  can  be  omitted. 

2)  If  the  signal  is  weak  compared  to  the  reverberation,  the  dominant 
noise  term  is  the  last. 

3)  A  coherent  detector  cross-correlates  the  output  as  a  single  receiver 
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against  a  properly  delayed  replica  of  the  transmitted  signal.  Its 
output  therefore  consists  of  the  first  two  terms  of  Equation  (l4) 
or  (15).  It  follows  that  the  instrumentation  of  Figure  1  behaves 
essentially  as  a  coherent  detector  when  the  signal  to  noise  ratio 
is  high.^-  When  the  signal  to  noise  ratio  is  low,  its  performance 
is  quite  different  and  will  be  seen  later  to  approximate  more  nearly 
that  of  an  incoherent  (power)  detector. 

4)  The  intermodulation  products  of  signal  and  reverberation  are  ef¬ 
fectively  range  gated.  In  Equation  ( lU )  only  scatterers  located  at 

c 

a  range  diflering  from  that  of  the  target  by  no  more  than  about 

c 

contribute  to  the  reverberhdorr  For  =  2nx50  (about  100  cps 

transmitted  bandwidth)  this  corresponds  to  a  range  variation  of 

about  16  feet.  In  Equation  (15)  scatterers  may  deviate  in  range 

1 

from  the  target  by  roughly  ^  agc.  With  a  transmitted  pulse  of 
0.5  sec  duration,  this  corresponds  to  a  range  variation  of  about 
1U5  feet. 

5)  A  somewhat  more  complicated  range  gating  effect  is  present  in  the 

terms  resulting  from  intermodulation  of  reverberation  with  reverber¬ 
ation.  According  to  Equation  (l4)  significant  contributions  to 
the  last  tem  ariBe  only  when  the  following  two  inequalities  sure 
satisfied . 

^It  will  be  shown  later  that  the  cross-correlation  between  the  second 
and  third  terms  of  Equations  (14)  and  (15)  is  small  under  most  conditions  of 
practical  interest.  Hence,  under  conditions  of  high  signal  to  noise  ratio, 
the  instrumentation  of  Figure  1  may  be  replaced  conceptually  by  a  coherent  de¬ 
tector  operating  in  a  background  of  reververation  with  twice  the  actual  power. 

p 

Note  that  t^  is  round  trip  time  to  scatterer  i.  Hence  a  range  increment 
of  s  feet  changeslt^  by  2s.  seconds . 
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(17) 


l(Vta)  -  <Tj-*b>l  <zr 

KW  +  (Tj'tb)l  <  2\(f2+  (18) 

The  spatial  region  to  which  these  relations  restrict  the  location  of 
contributing  scattered  pairs  is  indicated  by  the  crosshatched  region  of 


Figure  2.  Note  that  the  constraint  is  one  on  range  only.  Scatterer  pairs 
having  the  proper  range  relation  contribute  reverberation  regardless  of  their 
relative  bearing.'*’  The  origin  of  Figure  2  denotes  a  scatteirpair  such  that 
8catte»i  is  located  at  the  same  range  as  the  target  relative  to  receiver  a 
while  scatterer  J  is  located  at  the  saine  range  as  the  target  relative  to 
receiver  b.  The  total  range  variation  allowed  by  the  crosshatched  area  is 


If  they  are  simultaneously  illuminated  by  the  transmitting  beam. 

The  beam  parameters  are  introduced  later  by  assigning  angular  dependence  to 
a^  and  A. 
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clearly  of  the  order  of  V  2  c  \Ja^+a^/2  .  However  any  given  intermodulating 

pair  can  have  a  range  differential  no  greater  than  c/Ko  . 

c 

In  the  absence  of  frequency  modulation  [Equation  (15)3,  the  equivalent 
of  Equations  (17)  and  (18)  is 

iVa1  ’  (Tj*Vl  <V§oT  {19) 

and  |(ti-ta)  +  (Tj-tb)|  <  2  2  (20) 

The  region  of  contributing  scatterer  pairs  is  shown  in  Figure  3. 


6)  An  interesting  special  case  of  5)  [intermodulation  of  reverberation 
with  reverberation]  arises  when  i=J.  This  means  that  the  returns  from  a 
particular  scatterer  to  receivers  a  and  b  are  intermodulating  with  each  other. 
If  a  total  of  Ji  scatterers  is  illuminated,  the  number  of  terms  of  this  form 
is  N  whereas  the  number  of  modulation  products  of  different  scatterers  is 

p 

N  -N.  For  large  N  the  latter  predominate  and  the  former  are  of  practical 
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interest  only  if  the  return  from  one  or  more  scatterers  is  particularly  large. 
In  other  words,  the  case  i*J  is  essentially  that  of  false  targets.  For  i*j 
Equation  (17)  can  now  he  decomposed  into  the  tvo  expressions 


‘i  ■ T!  *  sr  *  <VV 

c 


Ti  - *i  *  s-  -  (W 

c 


Similarly  Equation  (18)  can  be  rewritten  as  follows 


(21) 

(22) 


tj  ♦  Tj  <  (ta«b)  ♦  2^0p2MT2ft'  (23) 

\  ♦  Ti > <  vV  - 2  2/2'  <2M 

ct,  is  the  distance  from  transmitter  to  receiver  a  via  the  1th  scatterer 
and  cTj  is  the  distance  from  transmitter  to  receiver  b  via  scatterer  i. 

It  follows  that  Equations  (21)  and  (22)  define  a  pair  of  hyperbolas  and 
Equations  (23)  and  (24)  a  pair  of  ellipses  between  which  false  targets 
must  lie  if  they  are  to  be  a  significant  source  of  confusion.  These 
constraints  on  false  target  locations  are  illustrated  in  Figure  4. 
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Figure  1*A  is  based  on  the  assumption  t  =t.  (target  broadside),  while 

&  D 

2 

Figure  UB  shows  the  case  ta~t^>  (target  significantly  clockwise  from 
broadside).  In  each  case  the  region  from  which  false  target  returns  may 
originate  (crosshatched  area)  is  a  pair  of  small  curvilinear  squares 
centered  at  the  target  location.  Thus,  as  far  as  false  targets  are  concerned, 
the  postulated  instrumentation  introduces  not  only  a  range  gating  effect 
but  also  a  gate  in  azimuth  and  elevation.*  In  practice  the  situation  is 
actually  likely  to  be  simpler  than  suggested  by  Figure  U,  because  the  target 
is  almost  certain  to  be  at  a  range  which  is  large  compared  with  the  spacing 
between  receivers  a  and  b.  In  that  case  the  ellipses  and  hyperbolas 
degenerate  into  circles  and  radial  straight  lines  respectively,  as  indicated 
in  Figure  5-  It  is  now  a  simple  matter  of  analytic  geometry  to  determine 

the  angular  spread  (Ag)  which  might 
contribute  a  false  target  return.  With 
a  true  target  broadside  (g*0)  one  obtains 


Figure  5 

l.e.  if  the  correlation  distance  of  the  transmitted  signal  in  the  water 
is  much  smaller  than  the  spacing  between  receivers,  Equation  (25)  reduces  to 

.  c/(Ka  ) 

Ag  ■  2 - j-2-  radians  (26) 

With  c/K°e  *  16  feet  (as  postulated  before)  and  d  *  160  feet,  this  leads 
to  an  effective  aperture  ag  *  0.2  radians.  Another  way  of  stating  this 

*The  significant  region  of  false  target  locations  in  3  dimensional 
space  is,  of  course,  the  solid  of  revolution  generated  by  rotating  the 
crosshatched  areas  about  the  axis  (a,b). 
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result  is  that  the  postulated  instrumentation  is  capable  (regardless  of  beam 

pattern  considerations)  of  resolving  two  targets  separated  in  angle  by 

about  0.1  radians.  An  expression  equivalent  to  Equation  (23)  for  f#0 

assumes  the  form  ./c/Ko  i  1/e/Ko  \ 

AB  *  sin  —  -  -  +  sin  B U  sin  -  -  sin  B  (27) 

where  each  arc  sine  is  interpreted  as  ir/2  when  its  argument  exceeds  unity. 

II.  Probability  Density  of  Scatterer  Location 

The  detector  output  z(t  +t  )  as  described  by  Equations  (14)  or  (15) 

a  D 

is  a  random  variable  whose  properties  are  dependent  on  scatterer  location 

to  the  extent  that  scatterer  position  influences  t^  and  Tj.  In  order  to 

describe  the  statistical  properties  of  z(t  +t.  )  it  is  therefore  necessary 

a  o 

to  postulate  a  statistical  model  of  scatterer  distribution  and  to  express 
the  properties  of  t^  and  Tj  in  terms  of  appropriate  paramenteia of  this  model. 
There  sure  two  major  sources  of  reverberation 

1)  Scatterering  centers  distributed  throughout  the  volume  of  water 
traversed  by  the  signal  (volume  reverberation). 

2)  Scatterering  centers  near  the  ocean  surface  or  bottom  (surface  or 
bottom  reverberation). 

For  the  purposes  of  this  report  the  simplest  possible  assumptions  vill 
be  made  concerning  the  distribution  of  each  type  of  scatterer. 

l)  Volume  reverberation.  Consider  a  large  volume  V  of  ocean 
surrounding  the  source  and  Including  all  of(but  encompassing  much  more  than) 
the  volume  illuminated  by  the  source.  The  location  of  the  ith  scatterer  will 
be  taken  as  a  random  variable  uniformly  distributed  over  V  and  independent 
of  the  location  of  any  other  scatterer.  Thus  the  number  N  of  scatterers  in 
V  has  a  Poisson  distribution.  Each  scatterer  is  assumed  to  be  stationary 
while  illuminated  by  the  transmitted  pulse.  The  scatterer  distribution 
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2 


is  described  most  simply  in  terms  of  a 
rectangular  coordinate  system.  Consider 
such  a  system  set  up  with  origin  at 
receiver  a  and  receiver  b  located  on 
■the  y  axis  as  shown  in  Figure  6.  The 
probability  that  the  ith  scatterer  lies 
in  [(x  ,y  ,z  ),(x.-t-dx,y.+dy,2.+dz)] 


Figure  6  is  clearly  (1/^dxdydz) .  Transforming 

this  result  into  a  spherical  coordinate  system  centered  at  receiver  a  as 
indicated  in  Figure  6,  one  obtains 


Pr  {ith  scatterer  in  (ri,ei,$^>(ri+dr,ej+de,$i+d$ )}  =  ^ri2sin©idrdQd|, 


28) 


Hence  the  probability  density  of  the  ith  scatterer  location  in  spherical 
coordinates  is 


P(r±  ,©i4»i )  = 


r.2sin©. 
V  1  l 


in  V 


(29) 


0  elsewhere 

If,  as  a  matter  of  computational  convenience  (and  without  influencing  the 
results  in  any  significant  way),  one  assumes  that  the  transmitter  is  located 
close  to  receiver  a  ,then 


t.  =  —  r. 
1  c  1 


(30) 


Hence  the  Joint  probability  density  of  t^,©^  and  ^  is 

lt  .  4  .  )  Bv  ‘i2  si"  ei  in  * 

p(ti,©i,4ii)  =  < 

0  elsewhere 

where  x  is  the  space  equivalent  of  V  under  the  transformation  (30). 

2)  Surface  reverberation.  The  scatterers  are  now  distributed  in  the 


(31) 
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immediate  vicinity  of  a  plane  (the  average  surface  or  bottom)  which  will 
be  taken  as  parallel  to  the  xy  plane  of  Figure  6.  If  Q  is  an  area  on  the 
surface  much  larger  than  that  illuminated  by  the  beam  and  if  one  assumes 
uniform  and  independent  distribution  of  scatterers  in  the  x  and  y  directions 
one  can  write 


Pr  {ith  scatterers  in  (x  ,y  ,z  ),(x  +dx  ,y  +dy  ,z.+dz  )}  =  i  exp 

1  1  1  111111  ^2ttA 

•  dx.dy.dz. 

l  i  l 


i  (vD)‘ 


2A 


(32) 


D  is  the  depth  of  the  array  relative  to  the  average  surface  and  A  is  the 
vertical  rms  spread  of  scatterer  locations.  The  Gaussian  distribution  in 
z  is  of  course  an  arbitrary  model.  Variations  in  the  parameter  A  can  be 
used  to  introduce  some  of  the  effects  of  surface  roughness.  Since  scatterer 
location  enters  the  expression  for  the  received  signal  only  via  the  time 
delays  t^  and  ,  it  is  clear  that  vertical  spreads  A  small  compared  with 
the  effective  illuminated  area  have  little  effect  on  the  return  signal. 

The  z  distribution  then  degenerates  effectively  into  a  delta  function  at 
z=D.  This  simplification  will  be  introduced  at  an  appropriate  point  later 
in  the  analysis .  It  should  also  be  pointed  out  that  the  primitive  model 
of  reverberation  used  in  this  report  ignores  the  time  dependence  of  scatterer 
locations . 

Transformation  into  spherical  coordinates  converts  Equation  (32)  into 
the  expression 


p(r. ,0. ,$. ) 

r  l  l  yi 


exp 


T  (r. cos0.-D)^1 

r 


r.  sine, 
l  i 


in  Q 


2A 

.i 

0  elsewhere 

which  is  the  equivalent  of  Equation  (29)  for  surface  reverberation.  Finally 
the  transformation  (30)  yields 


(33) 
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p(W>i)  *  i 


=3  1 
So  —  exP  - 
]l2vb 


!  7T  t.  cosS.  -D)^1 


2  i 


t . 2sin0 . 
1  l 


in  Y 


(3U) 


I.  2A 

0  elsewhere 

where  Y  is  the  space  equivalent  to  Q.  under  the  transformation  (30). 

The  Joint  distribution  of  t^,©^,  4k  is  sufficient  for  the  computation 
of  all  reverberation  statistics  because  the  random  variable  T.  can  be 

l 

expressed  in  terms  of  t^,©^,  and  4^  by  simple  geometrical  reasoning.  From 
Figure  6 


2,  2'  /,  2dyi“d 


Ri  *Vxi  +zi  +(Vd)  =7xi  +yi  +zi  I/1-  -r- 2 


xi  +yi  +zi 


2  ~  ri" 


dy.-dc/2 

r. 

i 


(35) 


where  use  has  been  made  of  the  assumption  d<<r^.  In  view  of  the  fact  that 
the  only  scatterers  contributing  to  the  reverberation  are  located  at  ap¬ 
proximately  the  same  range  as  the  target,  this  assumption  should  not  introduce 

any  significant  error. 

2 

Omitting  d  /( 2r as  negligible  compared  with  r^  and  transforming 

to  spherical  coordinates,  one  obtains 

R.  =  r.  +  d  sin©.  sin*.  (36) 

li  i  vi 

The  travel  time  of  sound  from  the  ith  scatterer  to  receivers  a  and  b 
is  r^/c  and  R^/c  respectively.  Thus,  assuming  the  source  to  be  located  near 


receiver  a 

T.  =  t,  +  —  sin0.  sin*, 
lie  l  ri 

III.  Average  Detector  Output 

The  DC  output  of  the  detector  is  the  expected  value  of  Equations  (1^) 
or  (15).  In  each  case  it  consists  of  a  signal  component 


(37) 


a  d 


(38) 
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and  a  noise  component.  The  second  and  third  terms  of  Equations  (ih)  and  (15) 

make  no  significant  contributions  to  the  average,  because  u  >>Ko>>-— 

e 

so  that  one  can  invoke  the  Riemann-Lebesgue  lemma.  For  the  same  reason 
all  terms  of  the  double  sum  for  which  i^j  can  be  ignored.  Thus  for  case 
A)  [bandwidth  determined  by  frequency  modulation] 


E{z(ta+tD»  Incise  "  E 

F  i 


-  2  ■'  2 


_i_  exp  <  j—  [(ti-ta)-(T.-t b)]‘ 


lL*i  Ti 


r 


exp 


[<ti-ta)*(T1-tb)r 

I 


’cos 


(  o  2  T 

[(ti-ta)-(Ti-tbm»o*f  \  WW1) 

^  F  (39^ 


T^  is  now  expressed  in  terms  of  t^,  ©^ ,  and  by  use  of  Equation  (37). 
Similarly  t^  can  be  written  as 

t  =  t  +  —  sin©^  sin<|>  (Uo) 

b  a  c  o  o 

where  and  ©q  measure  target  bearing  and  elevation  respectively.  For 

volume  reverberation  the  relevant  probability  distribution  is  given  by 

Equation  (31).  Substituting  Equations  (31),  (37),  and  (1*0)  into  Equation 

(39)  one  obtains  an  expression  which  must  be  integrated  over  the  variables 

and  ©j.  The  result  of  the  t^  integration  is  (after  considerable 

algebraic  manipulation) 
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The  second  exponential  in  Equation  (39)  differs  significantly  from  zero 
only  over  an  interval  roughly  equal  to  the  pulse  length, 
the  l^mitj  of  integration  can  be  taken  as  (-ao  ,ao)  and  a 
by  a  / t  without  introducing  appreciable  error. 


centered  at  2tg.  Hence 
^  /T^  can  be  replaced 
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Before  the  $  and  9^  integrations  can  be  performed,  it  is  necessary 

to  recall  that  the  coefficient  a^  determines  the  power  returned  from  the 
th 

i  scatterer.  is  therefore  dependent  on  the  beam  pattern  of  the  trans¬ 
mitter  and  receiver.  This  relationship  is  now  made  explicit  by  the  following 
expression 


bi  I  ( W  i 
sl^e.expr  TT 
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<V*o}  1 


(U2)1 


B_  and  B,  measure  the  pattern  width  in  the  9  and  direction  respectively 

til 

while  b^  is  the  scattering  cross-section  of  the  i  scatterer.  Equation  (42) 

2 

assumes  that  the  beam  pattern  is  centered  at  the  target. 

If  the  beam  dimer. sions  B,.  and  B.  are  small  compared  to  one  radian,  one 

9  9 

can  approximate  the  sinusoidal  functions  in  Equation  (Ul)  by  the  linear 
term  of  the  Taylor  Series  expansion. 

sin0.sin<|>.  *  sin9  sin4>  +  cos9  sin$  (9.-9  )  +  sin9  cos$  ($.-*)  .  (43) 

l  i  o  o  o  o  io  o  oio 

The  averaging  operation  in  extends  over  (0,it)  while  that  in  covers  a 

full  two  radians.  The  narrow  beam  patterns  already  postulated  enable  one 

to  extend  the  limits  of  integration  to  {-oo,  oo)  except  for  values  of  9q 

near  0  and  n  (straight  up  and  straight  down).  Since  these  particular 

directions  are  of  little  practical  interest,  the  complications  arising  near 

Ge=0  or  it  will  be  ignored, 
o 


1The  factor  ain9.  in  the  denominator  reflects  the  fact  that  the  $ 
dimension  of  a  fixed  element  of  area  varies  inversely  with  sin9.  Hence 
for  large  B  and  Equation  ( U2 )  corresponds  to  uniform  illumination  on 
any  sphere  about  tne  transmitter. 

2 

If  one  thinks  of  a  transmitter  operating  under  an  average  power 
constraint,  it  might  be  reasonable  to  postulate  that  the  total  radiated 
power  is  fixed  so  that  the  target  powtr  would  vary  inversely  as  B^B  . 

Such  an  assumption  is  easily  incorporated  into  the  analysis  by  dividing  the 
right  side  of  Equation  (42)  by  B^B  .  As  it  stands  Equation  (42)  implies  that 
the  power  returned  from  the  target  is  independent  of  Bq  and  B^. 
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The  actual  integrations  are  tedious  but  quite  straightforward.  The 
result,  after  extensive  algebraic  simplification,  is 


oo  7 

Idt  ij 

T  IT 

d9.  [dip.  z(t  +t_) 
i  i  a  D 

/  2 

!  p(t.,9i,$.)  =yir(cJF  +  -y-) 

o  3  i 

c  c 

2o„  W\ 

1  N  2v 

-oo  c 

)  -TT 

noise 

r 

\  i=r& 

exp<  - 


.2  2 

a  a) 

o 


°T  d  o 


Ue 

1  + 


t^2  2.2  2 

roTdp 


2  2  2  2  2  2  2 
where  p  '  =  B .  sin  9  cos  d>  +  Brt  cos  9  sin  $  • 

$  o  Yo  9  o  ro 

Equation  (hk)  contains  two  randan  variables,  the  scatterer  cross  section  b. 


and  N,  the  number  of  scatterers  in  the  volume  V.  Designating  the  average 
scatterer  cross  section  as  <b2>  and  postulating  an  average  of  K.  scatterers 
per  unit  volume  the  expected  value  of  the  detector  output  noise  becomes  finally 
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If  (K  Oj,  d  p  )/(8c  )  >>  1  the  argument  of  the  exponential  function 


becomes - - — 0  -  -2  § - rrrr  .^. — .  .  .  The  magnitude  of  this  ratio 

^2  2  transmitted  bandwidth 


is  almost  certain  to  be  large  compared  to  unity  so  that  the  DC  output  due 

p  2  2  2 

to  noise  is  negligible.  If  (Ko,pd  p  )/(8c  )  <<  1  the  argument  of  the 
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exponential  becomes  - 


For  ci=100  ft. 


=  2*  x  3500  rad/sec ,  and 


.2  2  2 
d  “o  P 

4c2 


»  u. 


2  2  2  2 

c  *  5000  ft/sec,  d  u>0  Ac  =  (70*)  .  Hence  the  exponential  is  of  an  order 
-10  2 

no  larger  than  e  for  o  i  0.0002.  Small  values  of  p  occur  near  0  =0, 

o 

$o*0,  and  0Q=  ♦  ♦_  &  .  Only  the  latter  combination  is  of  practical 

interest.  Postulating  B^«BQ=0.1  and  working  with  the  extreme  case,  0q=  —  , 
one  concludes  that  the  average  noise  output  is  very  small  unless  the  system 
is  trained  on  a  target  no  more  than  8°  from  the  endfire  direction.  The 
appearance  cf  Q  DC  component  of  noise  near  the  endfire  direction  is  reasonable, 
for  in  that  condition  the  difference  in  signal  travel  tine  to  the  two  receivers 
is  very  nearly  the  same  for  any  point  in  the  illuminated  volume.  In  other 
words,  because  of  poor  directivity  near  endfire,  the  detector  is  unable  to 
distinguish  between  a  point  target  and  a  volumetric  distribution  of  scatterers. 
In  effect  it  reports  the  accumulated  return  from  all  those  scatterers  which 
it  is  unable  to  resolve  as  on  equivalent  target.1 

The  derivation  of  the  average  noise  output  in  the  absence  of  frequency 
modulation  (case  B)  follows  the  same  pattern  and  leads  to  the  result 
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When  (d  q  )/(2c  'a^  )  <<  1  the  exponent  reduces  to  -(d  p  )/(4c  )  as 


Caution  is  in  order  in  any  quantitative  use  of  Equation  (46)  very  close 

to  the  endfire  direction,  for  the  use  of  only  linear  terms  in  Equation  (43) 

becomes  inappropriate  for  0  and  <4  within  Bx  and  B  of  (6  =  —  ,  A  =  x) . 

o  To  $  $  o  2  o  2 

Higher  order  terms  must  be  included  if  quantitatively  accurate  results  are 

desired. 
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under  the  equivalent  approximation  on  Equation  (46)  so  that  the  previous 
comments  concerning  DC  noise  output  near  the  endfire  condition  remain  ap¬ 
plicable. 

Equivalent  results  for  surface  reverberation  could  be  obtained  by  using 
the  distribution  of  Equation  (34)  in  place  of  Equation  (31).  The  computation 
is  quite  tedious  and  was  Judged  to  be  of  limited  interest  in  view  of  the 
anticipated  negligible  value  in  all  but  the  extreme  endfire  direction. 

IV.  Detector  Output  Variance 

In  order  to  evaluate  the  performance  of  the  detector,  it  is  necessary 
to  compute  the  output  fluctuation  as  well  as  the  average  output.  Since  the 
signal  component  of  Equations  (l4)  and  (15)  is  a  constant,  the  fluctuation  is 
entirely  due  to  the  noise.  For  the  frequency  modulated  case  [case  (A)]  one 
obtains  from  Equation  (lU) 
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Because  of  the  odd  number  of  indices  involved  in  terms©  and (6)  one 
can  immediately  invoke  the  Riemann-Lobesgue  lemma  and  assert  that  these 
terms  make  no  significant  contribution  to  the  mean  square  value.  Similar 
considerations  lead  to  the  conclusion  that  contributions  from  other  terms 
arise  only  for  the  fdllowing  combinations  of  indices, 
term©  :  j*J.' 
term©  :  i=i* 
term (3)  :  i=i '  ,J=J '  ,i£} 
i=J ,i'=J '  ,i^i' 
i-J’.J-i'.i* 


i=i'*  j=J ' 
term  ©  :  i=J 

The  actual  performance  of  the  averaging  operations  is  exceedingly 
cumbersome  but  basically  straightforward.  With  approximations  similar  to 
those  used  in  Section  III  one  obtains  in  the  case  of  volume  reverberation 
the  following  result  for  the  variance  of  z(t^+tD) 
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'Ifce  exponential  terms  art  all  similar  in  form  to  Equation  ( U6)  and 
are  therefore  negligible  except  very  near  to  horizontal  and  endfire.  When 
they  are  negligible  Equation  (1+9)  reduces  to 
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The  first  term  of  Equation  (50)  represents  the  power  of  the  intermodulation 
products  of  signal  with  the  returns  from  various  scatterers  [terms  and  (2) 
of  Equation  (1*8)].  The  second  term  results  from  the  intermodulation  of  the 
returns  from  different  scatterers  [term(D  i*i '  >Ja<)  *  »i?*J  in  Equation  (1+8)] 
while  the  last  term  gives  the  iotal  power  of  intermodulation  products  from  the 
same  scatterer  (false  target  power)  [term(|i  i*J*i,ssJ'  in  Equation  (1+8)]. 

The  radical  in  the  denominator  of  the  last  term  reflects  the  bearing  dis¬ 


crimination  against  false  targets  discussed  qualitatively  in  Section  I. 
Comparison  of  the  last  two  terms  reveals  that  the  last  term  is  likely  to  be 
small  compared  with  the  second  unless  the  scatterer  density  in  the  illuminated 
volume  is  very  low  (a  situation  in  which  the  reverberation  is  composed  of 
returns  from  a  few  strong  scatterers  -  which  may  clearly  be  regarded  as  false 
targets) . 

In  the  absence  of  frequency  modulation  a  completely  analogous 
computation  leads  to  the  follvoing  equivalent  of  Equation  (50) 
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When  the  primary  noise  source  is  surface  reverberation  p(t^,6^,<^)  is  given 
by  Equation  (3^).  The  averaging  operation  of  Equation  (U8)  with  respect 
to  this  distribution  is  basically  straightforward  but  even  more  cumbersome 
than  the  corresponding  computation  for  volume  reverberation.  It  has  been 
carried  out  only  for  terms®  ,  (2)  ,  and®  with  i=i’,  i^J-and  under  the 

assumption  that  the  vertical  dispersion  A  of  3catterers  satisfies 
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Equation  (52)  demands  that  A/sin©o  subtend  an  angle  small  compared  to 

Bq  at  the  source  [see  Figure  7],  The  approximation  implied  by  Equation  (52) 

& 


clearly  improves  with  range.  For  realistic  beam  widths  it  should  be  excellent 
except  at  extremely  small  ranges. 

With  the  stated  approximations  the  averaging  operations  lead  to  the 
following  result,  equivalent  to  the  first  two  (generally  dominant)  terms  of 
Equation  (50) 
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In  this  equation  Kg  is  the  number  of  surface  scatterers  per  unit  area, 
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Thus  tg  is  the  signal  round  trip  time  between  the  source  and  the  point 
on  the  surface  at  the  center  of  the  beam.  Similarly  is  the  differential 
in  round  trip  time  between  surface  returns  from  the  center  of  the  beam  and 
surface  returns  from  the  nominal  edge  (©o+Bq)  of  the  beam.  [See  Figure  7]. 

The  most  obvious  difference  between  surface  and  volume  reverberation 
as  expressed  by  Equations  (50)  and  (53)  is  the  presence  of  the  exponential 
terms  in  Equation  (53).  These  simply  reflect  the  fact  that  surface  rever¬ 
beration  if  a  problem  only  when  the  signal  round  trip  time  to  the  surface 
(nominally  t  )  is  close  to  the  signal  round  trip  time  to  the  target  (t  ). 


A  factor  of  the  form  exp 
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term  has  been  omitted  as  negligible  under  most  reasonable  operating  conditions, 
particularly  when  D^d. 
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It  is  interesting  to  observe  that  the  allowed  time  differential  t  -t  is  of 

as 

the  order  of  oQ  for  the  noise  component  resulting  from  intermodulation  of 

signal  with  reverberation  while  the  differential  t  -t  can  be  as  large  as 

as 

/  2  2  2' 

WOp  +o^  /2+Og  for  the  intermodulation  products  of  reverberation  with  rever¬ 
beration.  Thus  the  use  of  very  long  pulses  may  create  a  surface  reverberation 
problem  at  ranges  where  it  would  not  otherwise  exist  if  the  dominant  component 
of  Equation  (53)  for  t  =t  is  the  second  term  (low  signal  to  noise  ratio). 

0.  S 

This  difficulty  does  not  arise  when  the  first  term  dominates  (high  signal  to 
noise  ratio).  The  physical  explanation  for  this  phenomenon  is  simple :He  frequenc- 
modulation  imposes  a  narrow  range  gate  on  reverberation  (see  Section  I). 

Only  scatterers  located  such  that  the  travel  time  of  the  signal  to  them  is 
almost  exactly  the  same  as  the  travel  time  to  the  target  can  contribute  to 
the  intermodulation  product  of  signal  and  reverberation.  Whether  the  group 
of  soatterers  at  suitable  ranges  does,  in  fact,  contribute  reverberation 
depends  only  on  whether  they  lie  within  the  illuminated  surface  area.  This 
is  purely  a  matter  of  geometry  and  independent  of  the  pulse  length.  Inter- 
modulation  of  reverberation  with  reverberation,  on  the  other  hand,  occurs 
for  any  scatterer  pair  separated  from  each  other  by  no  more  than  the  effective 
range  gate.  Whether  a  particular  combination  contributes  reverberation  at 
the  observation  time  (t&+tp)  therefore  depends  both  on  whether  it  received 
illumination  and  on  when  it  received  illumination,  i.e.  on  both  geometry  and 
pulse  length. 

Another  interesting  point  of  comparison  between  volume  and  surface 

reverberation  is  the  radical  in  the  second  terms  of  Equations  (5r)  and  (53) 

respectively.  When  o  >>  (a„  +o~/2)  the  radical  in  Equation  (53)  reduces 

V  r  1 
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to  *^»Op  *Oip  /2.  Thus  the  reverberation  power  is  limited  by  the  pulse 

2  2  2 

duration,  as  in  the  case  of  volume  reverberation.  When  <j  <<(av  +om  /2) 

y  FT 

the  radical  in  Equation  (53)  reduces  to  V"2 eg.  Thus  further  increases  in 
pulse  length  do  not  increase  the  peak  surface  reverberation  power.  This  is, 
of  course,  reasonable  because  the  pulse  duration  is  now  sufficiently  large 
for  reverberation  from  the  entire  illuminated  surface  area  to  reach  the  receive 
simulaneously. 

Additional  differences  between  volume  and  surface  reverberation  concern 

the  t  dependence  (t  and  t  ”  with  volume  reverberation  versus  t  ~  and 
a  r  a  a  a 

tQ-^  with  surface,  reverberation)  and  the  appearance  of  the  factors  sin9Q 
2 

and  sin  0  in  the  case  of  surface  reverberation.  These  differences  are 
o 

attributable  to  the  different  rate  of  growth  with  range  of  the  number  of 
scatterers  contributing  reverberation  in  the  two  cases. 

V.  Output  Signal  to  Nbise  Ratio 

The  effectiveness  of  the  detector  is  characterized  by  the  relation  of 
its  output  signal  to  noise  ratio  with  the  signal  to  noise  ratio  existing  in 
the  water.  The  output  signal  power  is  simply  the  square  of  the  first  term 
of  Equation  (13).  Using  t  *t.  one  obtains  from  Equation  (13)  and  the  first 
two  terms  of  Equation  (50)  (volume  reverberation,  frequency  modulated  signal) 
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A  simple  computation  from  Equation  (U)  gives  the  average  reverberation 


*The  effective  integration  time  o?  of  the  low-pass  filter  should  clearly 
be  of  the  order  of  the  pulse  width  o_.  A  more  precise  optimum  will  be  obtained 
later. 
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(57) 


power  at  eaob  receiver. 

Average  reverberation  power  »(j)2  g— o^, 

t 

a 

If  t^t^  the  peak  signal  power  at  each  receiver  is  t?/t\ 


Hence 
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Thus  Equation  (56)  con  be  rewritten  in  the  form 
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Several  features  of  this  expression  are  interesting. 
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Using  the  definition  of  oc  [Equation  (10)]  it  is  a  simple  matter  to  determine 
the  optimum  <jp  for  either  case.  From  Equation  (60)  one  obtains 


T 

°F 


(62) 


The  optimum  in  Equation  (6l)  is  infinite.  However,  it  is  clear  from 
r 

Equation  (10)  that  only  small  gains  in  output  signal  to  noise  ratio  can  be 

2  2 

made  by  increasing  Op  above  <£  /2.  Equation  (62)  therefore  gives  a  suitable 
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value  of  o for  all  Input  signal  to  noise  ratios.  With  this  choice 


2  *  2,, 

Op  +oT  / 2 


(63) 


so  that  Equations  (60)  and  (61)  nay  now  be  interpreted  as  follows: 

When  the  input  signal  to  noise  ratio  is  less  than  unity,  the  output 
signal  to  noise  ratio  varies  with  the  square,  of  the  input  signal  to  noise 
ratio,  a  type  of  behavior  generally  associated  with  incoherent  detection. 

When  the  input  signal  to  noise  ratio  exceeds  unity,  the  output  signal  to  noise 
ratio  varies  linearly  with  the  input  signal  to  noise  ratio,  a  characteristic 
typical  of  coherent  detection  schemes.  Equation  (61)  is,  in  fact,  simply 
one-half  of  the  output  signal  to  noise  ratio  of  a  coherent  detector  (de¬ 
modulating  the  output  of  &  single  receiver  against  a  delayed  replica  of  the 
transmitted  signal)  .  The  factor  of  one-half  is  due  to  the  presence  of  noise 
in  each  channel  of  the  correlation  detector.  The  replica  used  in  the  coherent 
detector  is,  of  course,  noise-free.  Equation  (59)  with  Op“oT/2,  Kot"2it  x  100, 
and  ot»0.5  is  plotted  in  Figure  (8).  Also  shown  is  the  corresponding  curve 
for  the  coherent  detector.  If  one  chooses  an  output  signal  to  noise  ratio 
of  6  db  as  the  mlmimum  level  at  which  some  significant  detection  capability 
may  be  said  to  exist,  Figure  (8)  indicates  that  the  correlation  detector  requires 
about  8  db  more  input  signal  to  noise  ratio  than  the  coherent  detector  to  achieve 
this  minimal  performance.  Even  if  one  penalizes  the  coherent  detector  by  3 
db  on  the  grounds  that  the  postulated  knowledge  concerning  carrier  phase  is 
almost  certainly  unavailable  in  practice,  there  remains  a  differential  of 
about  5  db. 


The  presence  of  ambient  noise  -  ignored  thus  far  -  would  also  mitigate 
agaitPt  a  choice  of  effective  smoothing  times  Op  much  larger  than  the  duration 
of  the  demodulated  signal  pulse. 
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The  factor  o^,(Koc)  appearing  in  Equation  (59)  can  be  rewritten 

W  *  =  I  ITT^T  ^ulth  °f  =  V^)-  < 6k) 

Since  Ko^,  is  a  direct  measure  of  the  signal  bandwidth,  its  inverse  may¬ 
be  interpreted  as  the  correlation  time  of  the  signal  (and  hence  of  the 
reverberation).  The"processing  gain"  factor  ( 6U )  is  therefore  the  ratio  of 
signal  duration  to  signal  correlation  time.  It  can  be  made  substantial  by 
the  use  of  wide  signal  bandwidths.  It  is,  of  course,  nothing  more  than  the 
quantitative  expression  of  the  range  gating  effect  discussed  earlier. 

In  the  absence  of  frequency  modulation  the  equivalent  of  Equation  (59)  is 
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Vo 
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T 

i/T 


(65) 


(66) 


the  "processing  gain"  is  of  the  order  of  unity.  All  other  comments  concerning 
Equation  (59)  apply  to  Equation  (65)  also. 

In  order  to  avoid  misinterpretation  of  Equations  (59)  and  (65)  it  should, 
perhaps,  be  pointed  .out  that  the  input  signal  to  noise  is  also  a  function  of 
o^.  Thus,  from  Equation  (58) 


(-) 

Vin 


JL 

o„ 


(67) 


Hence  in  the  presence  of  frequency  modulation  [Equation  (59)) 
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Tn  the  absence  of  frequency  modulation  [Equation  (65)] 
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The  use  of  short  pulses  is  advantageous  in  the  absence  of  frequency  modulation 
because  the  bandwidth  (and  hence  the  range  gating  effect)  is  determined 
exlusively  by  oT-  In  the  presence  of  frequency  modulation  pulse  length 
is  a  secondary  factor,  for  the  rise  in  bandwidth  offsets  the  decrease  in 
input  signal  to  noise  ratio  accompanying  an  increase  in  oT>  For  fixed  band¬ 
width  (fixed  Ka,p)  the  output  signal  to  noise  ratio  increases  monotonically 
as  decreases.  Thus  the  use  of  short  pulses  (with  large  values  of  K  to 
maintain  constant  bandwidth)  is  indicated.  However,  once  the  input  signal 
to  noise  ratio  exceeds  unity,  further  decreases  in  oT  cause  only  a  small 
increase  in  output  signal  to  noise  ratio. 

In  the  case  of  surface  reverberation  the  signal  to  noise  ratio  is  critically 

2  2 

dependent  on  the  exponential  factors  in  Equation  (53).  When  (t  -t  )  <<  oQ 

a  s  w 

both  of  these  factors  are  close  to  unity  and  one  obtains  the  following 
equivalent  of  Equation  (59)  (frequency  modulated  signal) 
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With  Op  *  this  result  reduces  to 
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When  oT  «  0g 
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so  that  the  "processing  gain"  is  the  same  as  with  volume  reverberation.  On 
the  other  hand,  when  o^,  » 


oT(Kae) 


a9^T2+ae2‘y 
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Thus  the  "processing  gain"  rises  with  increasing  pulse  length  to  a  maximum 
determined  by  the  geometry  of  the  illuminated  surface  area.  As  in  the  case  \ 
of  volume  reverberation,  this  <5^  dependence  is  somewhat  deceptive  because 
of  the  <JT  dependence  of  the  input  signal  to  noise  ratio.  For  surface  rever¬ 
beration 
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Thus  the  output  signal  to  noise  ratio  increases  with  pulse  duration  only 
when  o_  exceeds  o0,  i.e.  the  entire  surface  area  within  the  beam  receives 
illumination  simultaneously. 

s  s 

For  high  input  signal  to  noise  ratios  [(jj')Q  “  (j^)in]  one 
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for  all  values  of  oT/o0.  These  observations  are,  of  course,  also  evident  from 
Equation  (53) . 

When  a_  »  oQ  and  the  target  is  so  located  that  a.  «  t  -t  « 

T  B  das 


2  2  21 

°F  +0T  ^+aB  t*ie  ^*rst  exponential  term  in  Equation  (53)  vanishes  while  the 

the  second  assumes  a  value  close  to  unity. 

In  that  case 


fix  „  I  ewT 
Vo  4 


'i 


l/a_  /2+o_  +o 


T 


6 


|/(aT2/2+oF2)'(oT2/2+aQ2) 


(K°c>  <f>in 


or  with  oF  -  ot//7 


/S  x  .  1  ,„n  x  bVt  +08  /Sv  2 

(No}  4  °T(Koc)  2.  2  N  in 

o  /2+o 


(77) 


(78) 


Thus,  it  is  possible  to  find  operating  conditions  under  which  the  output  signal 
to  noise  ratio  varies  with  the  square  of  the  input  signal  to  noise  ratio  for 
all  values  of  input  signal  to  noise  ratio. 

VI.  Comparison  With  Detector  Operating  in  Ambient  Noise  Limited  Environment 
Throughout  t..e  preceding  discussion  the  received  noise  was  assumed  to 
consist  exlusively  of  reverberation.  It  will  be  interesting,  for  purposes  of 
comparison,  to  evaluate  now  the  performance  of  the  same  receiver  in  a  noise 
environment  consisting  predominantly  of  ambient  noise.  The  ambient  noise  will 
be  regarded  as  a  stationary  Gaussian  random  process  with  the  power  spectrum. 
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or  the  autocorrelation  function 
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One  may  think  of  Nq  as  the  spectsl  level  of  a  broadband  noise  whose 
spectrum  is  being  shaped  by  a  bandpass  filter  of  halfwidth  centered  at 
frequency  In  practice,  this  filter  might  be  the  inherent  frequency  response 

characteristic  of  the  receiving  hydrophones  or  it  might  be  a  narrower  filter 
inserted  into  the  hydropoone  outputs  to  discriminate  against  frequencies  outside 
of  the  signal  range. ^ 

If  one  designates  the  noise  appearing  at  the  two  receiver  outputs  as 
na(t)  and  n^(t)  respectively,  then  from  Figure  1  and  Equation  (1) 
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The  detector  output  z(t)  is  related  to  y(t)  through  Equation  (6), 
the  low  pass  filter  weighting  function  w(t)  being  specified  by  Equation  (7). 
The  signal  component  of  z(t)  is  the  same  as  in  the  reverberation  limited  case. 
With  tr«ta~tb  the  noise  component  assumes  the  form 
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xAny  such  filter  is  assumed  to  be  sufficiently  broad  so  that  it  will  have 
negligible  effect  on  the  signal. 
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If  receivers  £  and  £  are  separated  by  more  than  a  very  few  feet,  it  is 
probably  reasonable  to  assume  that  nfl(t)  and  nb(t)  are  statically  independent. 
It  also  appears  reasonable  to  assume  that  each  noise  component  has  zero  mean. 
In  that  case  the  expected  value  of  z(t)  is  zero  and  its  variance  is  simply 
the  sum  of  the  mean  square  values  of  the  three  terms  in  Equation  (82) .  Thus 
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where  R(t)  is  given  by  Equation  (80)  .  The  first  two  terms  differ  only  in 
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the  attenuation  factor  A/t  vs  A/t, 

a  b 


Since  t  ■  t.  for  all  cases  of  interest 
a  b 


this  difference  will  be  ignored.  The  averaging  operation  can:  now  be  performed 
without  difficulty.  After  algebraic  simplification  the  result  is 
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The  output  signal  power  is  still  given  by  the  square  of  the  first  term 
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of  Equation  (13)  .  Hence 
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The  input  signal  to  noise  ratio  at  each  receiver  is  easily  computed  to  be 
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The  output  signal  to  noise  ratio  can  now  be  written  in  a  form  equivalent 
to  that  of  Equations  (59)  and  (65) 
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In  order  to  obtain  a  comparison  with  the  reverberation  limited  case,  it 
is  necessary  to  select  a  suitable  value  for  the  noise  spectral  width 
The  best  value  for  this  purpose  is  that  which  gives  the  ambient  noise  spectrum 
the  form  of  the  reverberation  spectrum. 

The  reverberation  autocorrelation  can  be  obtained  from  Equation  (3) 
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Carrying  out  the  averaging  operation  over  all  random  variables  one  obtains 
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Thus  the  reverberation  autocorrelation  function  ^Equation  (89)j  and  the 
ambient  noise  autocorrelation  function  are  identical  in  form  if 
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It  is  a  simple  matter  to  demonstrate  that  Equation  (91)  reduces  to  Equation 
(59)  when  the  bandwidth  of  the  transmitted  signal  is  determined  primarily 
by  the  frequency  modulation  Jcase  (A)j  .  In  the  absence  of  frequency  modulation 
|^case  (B)j  Equation  (91)  reduces  to  Equation  (65).  There  is  no  reason  to  doubt 
that  a  similar  check  could  be  obtained  for  frequency  modulated  signals  of  a 
bandwidth  comparable  to  l/0^,jfor  which  no  formal  computations  were  carried 
out  in  the  reverberation  limited  case.  Thus  it  appears  that  Equation  (91) 
may  be  regarded  as  a  generalization  of  Equations  (59)  and  (65) ,  covering 
intermediate  values  of  frequency  deviation  KoT  as  well  as  the  extremes  represent 
ed  by  cases  (A)  and  (B) .  One  therefore  concludes  that  the  detector  responds 
in  identical  fashion  to  reverberation  and  to  ambient  noise  of  the  same  power 


If  is  indeed  the  halfwidth  of  a  bandpass  filter  as  suggested  earlier, 
then  choice  of  the  value  given  by  Equation  (90)  would,  of  course,  result  in  some 
modification  of  the  signal  spectrum. 
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and  spectral  shape  as  long  as  the  assumptions  common  to  Equations  (59)  and 
(65)  are  valid.  These  are  basically  the  assumptions  permitting  the  reduction 
of  Equation  (49)  to  the  first  two  terms  of  Equation  (50).  In  physical  terms 
they  are 

1)  No  significant  percentage  of  the  total  reverberation  power  must 
be  contributed  by  any  small  collection  of  scatterers  with  large 
cross-section  (false  targets) . 

2)  The  parameter  p  a  combination  of  looking  angle  and  beam  dimensions 

9 

given  by  Equation  (45),  must  be  sufficiently  large  so  that  accumulationr 
of  unresolvable  small  scatterers  do  not  contribute  a  significant  per¬ 
centage  of  the  total  reverberation  power.  Practically  speaking, 
this  is  likely  to  be  a  problem  only  for  narrow  beams  directed  very 
close  to  horizontal  and  endfire. 

Requirements  1)  and  2)  are,  of  course,  precisely  the  conditions  under 
which  one  can  invoke  the  Central  Limit  Theorem  and  argue  that  the  combined 
returns  from  many  small  scatterers  form  a  Gaussian  random  process  whose 
Statistical  properties  are  then  completely  determined  by  its  autocorrelation 
function.*' 

The  above  remarks  should  not  be  interpreted  to  mean  that  there  is  no 
significant  difference  between  reverberation  and  ambient  noise,  even  when 
conditions  1)  and  2)  are  satisfied.  What  they  do  mean  is  that  any  differences 
in  performance  are  due  to  differences  in  the  input  signal  to  noise  ratio. 

Perhaps  the  most  important  feature  of  the  input  signal  to  noise  ratio  in  the 
ambient  noise  case  is  its  inverse  dependence  on  0^  j^see  Equation  (86) j  ,  If 

*See,  for  instance,  H.L.  Van  Trees, "Optimum  Signal  Design  and  Processing 
for  Reverberation  Limited  Environment'.'  IEEE  Transactions  on  Military  Electronic.' 
Vol.  MIL  9,  pp.  212-229,  July,  Oct.  1965. 
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Is  adjusted  in  accordance  with  Equation  (90),  i.e,  in  accordance  with 
the  signal  bandwidth,  one  does  not  obtain  the  Improvement  of  output  signal  to 
noise  ratio  with  increased  bandwidth  (frequency  modulation)  observed  in  the 
reverberation  limited  case.  Thus  Equation  (60)  ^low  input  signal  to  noise 
ratioj  now  varies  inversely  with  bandwidth  (Ko^)  while  Equation  (61)  j^high  input 
signal  to  noise  ratioj  is  independent  of  (Ko^,)  .  The  use  of  narrowband  trans¬ 
mitted  signals  is  therefore  indicated  for  ambient  noise  limited  operation, 
provided  that  Doppler  shift  considerations  (which  have  been  ignored  in  this 
analysis)  permit  the  use  of  a  sufficiently  narrow  filter  bandwidth 
VII.  Concluding  Remarks 

One  of  the  primary  purposes  of  this  report  is  to  provide  an  analytical 
framework  or  future  studies  of  various  problems  in  active  sonar  detection. 

Some  of  the  specific  assumptions  used  here  (stationary  target  and  scatterers, 
independence  of  location  of  different  scatterers,  absence  of  distortion  or 
multipath  effects  in  the  return  signal)  are  clearly  unrealistic  in  many  cases 
and  will  have  to  be  modified  in  subsequent  investigations.  In  the  meantime 
the  detection  problem  has  been  analized  at  least  under  highly  idealized  conditio 
and  insight  has  been  gained  into  such  questions  as  the  conditions  under  which 
reverberation  can  be  regarded  as  equivalent  to  Gaussian  ambient  noise  with  a 
known  spectrum.  The  results  of  the  present  analysis  should  therefore  permit 
simplifications  of  later  analyses  employing  more  realistic  models  of  the 
transmission  and  reverberation  mechanisms. 

The  first  steps  towards  the  generation  of  a  better  model  should  probably 
be  the  inclusion  of  Doppler  shifts  due  to  source,  target,  and  scatterer  motion 
and  recognition  of  the  fact  that  the  target  return  may  not  be  simply  a  delayed 
and  Doppler  shifted  replica  of  the  transmitted  signal.  Work  along  these  lines 
is  now  in  progress. 
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